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Preface 


Interest  in  the  use  of  solar  energy  to  supplement 
conventional  methods  of  heating  the  interior  of  structures 
has  grown  with  the  price  of  conventional  fuels.  A  method  of 
predicting  the  thermal  performance  of  structures  heated 
passively  by  the  sun  would  be  very  useful  in  determining 
the  cost  effectiveness  of  a  particular  design.  A  computer 
simulation  model  capable  of  estimating  the  hourly  perform¬ 
ance  of  a  residential  size  structure  with  a  variety  of 
passive  solar  heating  features  is  described  in  this  work. 
Many  designs  can  be  investigated  by  making  changes  in  the 
input  parameters. 

I  wish  to  thank  Dr.  James  1.  Hitchcock,  my  thesis 
advisor,  for  the  guidance  and  encouragement  he  gave  to  me. 
Without  his  efforts,  this  document  would  not  exist. 

I  would  especially  like  to  thank  my  wife  Helen.  Her 
untiring  support  and  understanding  during  this  busy  time 
helped  me  more  than  she  will  ever  know. 
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. — A  simple  lumped  capacitance-resistance  model  is  used 
to  simulate  heat  flow  in  a  residential  size  structure  heated 
passively  by  the  sun.  The  model  takes  the  form  of  an  analo¬ 
gous  electrical  circuit.  A  computer  program  was  written  to 
analyse  the  circuit.  By  altering  the  input  parameters  of 
the  program,  the  thermal  performance  of  a  wide  variety  of 
passive  solar  designs  can  be  investigated  for  any  geographi¬ 
cal  location.  By  comparing  program  generated  data  to  data 
taken  from  experimental  test  cells  in  Los  Alamos,  New  Mexico, 
it  was  found  that  the  simulation  program  predicted  energy 
use  to  within  4  percent  of  measured  values.  Also,  the 
computer  program  predicted  temperature  Bwings  to  within 
16  percent  of  measured  swings.  Correlation  with  empirical 
methods  of  calculating  monthly  and  annual  savings  in  fuel 
use  for  heating  was  poor.  Using  the  simulation  calculations 
as  a  base,  the  predictions  of  annual  savings  differed  by  as 
much  as  76  percent.  — 
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^  THE  SIMULATION  OF  A  PASSIVE  SOLAR  ENERGY  SYSTEM 

I  Introduction 

Background 

Only  recently  has  interest  in  solar  energy  for  space 
heating  become  widespread.  From  ancient  civilizations,  how¬ 
ever,  comes  evidence  of  construction  with  the  thermal  effect 
of  the  sun  in  mind.  Early  Greeks,  Egyptians,  and  Persians 
built  structures  that  accepted  solar  radiation  in  the  winter, 
yet  shaded  the  structure's  interior  in  the  summer.  In  North 
America,  the  Pueblo  Indians  built  cliff  dwellings  in  the 
Southwest  as  early  as  1000  A.D.  Erected  in  the  depressions 
f  of  vertical  south  facing  natural  formations,  the  buildings 

were  shaded  from  the  summer  sun  by  overhangs.  Wintertime 
brought  the  sun  to  lower  elevation  angles,  where  the  adobe 
walls  were  directly  irradiated.  These  bricks  held  their 
warmth  long  after  the  sun  went  down.  Whether  the  Pueblos 
first  erected  their  dwellings  for  thermal  comfort  or  as 
defense  from  enemies  is  uncertain.  Nonetheless,  the  bene¬ 
fits  of  passive  solar  control  were  certainly  recognized 
and  exploited  in  later  constructions. 

In  this  century,  solar  energy  has  been  largely  ignored. 
Vast  quantities  of  coal,  oil,  gas,  and  water  power  at  low 
cost  have  made  the  added  cost  of  passive  solar  construction 
unnecessary.  Now,  with  utility  costs  rising,  efficient 
solar  designs  can  be  more  economical  in  the  long  run  by 
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offsetting  higher  construction  costs  with  lower  utility  bills. 

Economic  analysis  methods  exist  to  determine  the  ad¬ 
visability  of  selecting  a  specific  solar  design*  given  its 
thermal  performance.  Calculating  thermal  performance  is 
more  difficult.  Around  the  country*  climatic  conditions! 
temperatures,  precipitation*  solar  insolation*  cloudiness, 
and  air  quality  vary  widely.  These  factors  and  more  affect 
how  well  a  solar  house  will  perform.  Procedures  to  calcu¬ 
late  thermal  performance  began  to  emerge  in  the  1960's. 

The  two  basic  types  of  procedures  are  empirical  methods  and 
simulation. 

Empirical  techniques,  such  as  the  popular  Solar 
Savings  Fraction  (SSF)  method  developed  by  Balcomb  (Ref  1), 
allow  one  to  estimate  the  savings  in  monthly  heating  bills 
for  a  specific  solar  design,  based  on  a  few  correlating 
parameters.  Use  of  the  SSF  method  entails  balancing  monthly 
solar  heat  gain  against  monthly  heat  loss  from  a  structure, 
then  reading  from  a  graph  the  percentage  that  furnace  use 
will  decrease  compared  to  a  similar  house  without  a  solar 
collector.  Calculations  take  only  a  few  hours  and  results 
are  accurate  for  long  term  analysis.  The  problem  with 
empirical  techniques  is  that  safeguards  against  poor  designs 
are  minimal.  For  example,  increasing  the  solar  collector 
area  may  very  well  cut  down  on  heating  bills,  but  may  also 
cause  the  interior  to  overheat  on  a  sunny  day. 

To  guard  against  hidden  design  defects,  one  must  look 
at  the  temperature  swing  (the  difference  between  the 
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maximum  and  minimum  interior  temperature)  for  the  structure. 
This  involves  a  great  deal  of  lengthy  simulation  calcula¬ 
tions  (solving  sets  of  differential  equations)  and  is  done 
on  the  computer.  Several  heat  transfer  simulation  programs, 
such  as  the  Transient  System  Simulation  Program  (TRNSYS ), 

(Ref  2),  and  D0E2  (Ref  3)  have  been  developed  recently. 

These  large,  complex  programs  calculate  heat  flow  rates 
for  both  active  and  passive  solar  heating  systems.  At  hourly 
or  shorter  intervals,  the  programs  will  report  temperatures 
in  various  building  components.  Using  the  print-out,  the 
designer  can  determine  not  only  the  economy  of  a  design, 
but  can  also  insure  that  interior  temperatures  stay  within 
acceptable  limits.  Such  programs  do  have  drawbacks,  however. 
While  versatile,  they  are  large.  The  1974  version  of  TRNSYS 
consists  of  almost  5000  FORTRAN  statements.  Later  versions 
are  larger.  In  addition,  these  large  programs  can  be 
cumbersome  to  operate,  especially  if  a  variety  of  designs 
are  under  consideration. 

Pu*3>g§e 

The  purpose  of  this  study  was  to  devise  an  efficient 
transient  simulation  program  capable  of  analysing  a  var¬ 
iety  of  passive  solar  designs.  It  should  give  the  energy 
savings  for  a  particular  design  and  predict  with  reasonable 
accuracy  the  temperature  swing.  The  program  should  be  small 
enough  to  fit  on  the  micro-computers  that  are  showing  up  in 
design  offices  and  be  quick  and  easy  to  use.  Combined  with 
an  eoonomic  analysis  method,  the  solar  program  could  be 
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used  to  pick  the  most  economical  design  from  several  such 
designs. 

Approach 

The  structure  was  modeled  using  an  electrical  circuit 
analog.  Resistors  and  capacitors  were  used  to  simulate 
insulation  and  thermal ’ mass .  respectively.  Electrical  po¬ 
tentials  corresponded  to  temperatures  at  physical  locations 
within  the  structure.  Differential  or  algebraic  equations 
for  each  node  in  the  circuit  diagram  were  derived  and  these 
equations  were  solved  on  the  computer  using  numerical  tech¬ 
niques.  Output  gave  a  history  of  solar  heat  gain,  nodal 
temperatures,  and  auxiliary  heating  requirements  for  a 
twenty  four  hour  period  on  typical  days  for  each  month  of 
the  year. 

The  development  of  the  transient  simulation  program 
is  presented  in  the  next  two  sections.  Then,  the  program 
verification,  performed  by  attempting  to  accurately  simulate 
an  instrumented  test  structure  set  up  in  Los  Alamos,  New 
Mexico,  is  presented.  Finally,  a  hypothetical  residential 
sise  structure  located  in  Columbus,  Ohio  was  modeled.  The 
thermal  performance  of  a  conventional  structure  was  compared 
to  the  performance  of  several  passive  solar  alternatives. 


II  Solar  Design  Analysis 


The  details  a  solar  engineer  must  consider  in  designing 
a  structure  that  will  he  heated  by  the  sun  are  presented  in 
this  section.  First,  the  design  concepts  involved  will  be 
discussed.  Next,  heat  flow  in  a  structure  will  be  examined. 
Thermal  energy  for  the  building  comes  from  the  sun,  from 
internal  equipment  and  people,  and  from  a  furnace.  Heat  is 
lost  through  the  ceiling,  walls,  and  floor.  Infiltration 
of  outside  air  in  and  inside  air  out  also  contributes  to 
heat  loss.  After  detailing  the  heat  flow,  the  Life  Cycle 
Savings  Method  of  economic  analysis  will  be  described. 

Finally,  a  few  human  factors  will  be  considered.  In  order 
to  design  a  comfortable  house,  one  must  establish  a  defini¬ 
tion  for  "comfortable." 

Solar  House  Designs 

Basically  two  categories  of  solar  heating  systems  exist. 
An  active  system  involves  the  use  of  pumps  and  or  fans  to 
move  a  working  fluid,  air  or  water,  from  the  solar  energy 
collector,  usually  mounted  on  the  roof,  to  a  separate  storage 
mass  for  heat  retention  or  directly  to  the  living  area.  With¬ 
out  some  form  of  energy,  such  as  electricity,  to  run  the  pumps 
or  fans,  an  active  solar  energy  system  will  not  work.  A  pas¬ 
sive  solar  heating  system  requires  no  auxiliary  energy  source. 
The  structure  itself  is  designed  to  collect  and  store  solar 
energy.  The  walls  and  floor  can  do  double  duty,  supporting'* 
the  roof  and  storing  heat.  Large  windows  not  only  allow  one 
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Pig.  1  Structure  with  Direct  Gain  Solar  Heating 

to  see  out,  but  admit  solar  radiation.  All  heat  flow  in  the 
passive  solar  structure  is  by  radiation,  conduction,  and 
natural  convection.  Passive  system  designs  can  be  further 
broken  down.  One  can  consider  direct  gain  or  indirect  gain 
configurations . 

The  simplest  passive  solar  heating  concept  is  direct 
gain  heating,  as  shown  in  Fig  1.  The  sun  shines  through 
south  facing  glazingi  provided  the  structure  is  located  in 
the  Northern  Hemisphere i  straight  into  the  living  area. 

The  floors,  walls,  and  furniture  heat  up  as  a  result  of 
this  direct  insolation.  The  interior  and  occupants  are 
warmed  by  convection  and  radiation  from  the  heated  surfaces. 
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The  living  area  acts  as  solar  collector ,  storage  mass,  and 
heat  distribution  system. 

'  The  advantages  of  a  direct  gain  system  are  its  simplicity, 
low  cost,  and  large  window  area.  As  shown  in  Fig  1,  the  sys¬ 
tem  consists  of  windows  and  thermal  mass.  These  are  easy 
and  relatively  inexpensive  to  build.  In  addition,  the 
windows  allow  pleasing  visual  access  to  the  outside. 

Direct  gain  systems  also  have  disadvantages.  The  large 
expanse  of  glass  can  allow  the  interior  to  become  too  bright 
on  a  sunny  day.  Also,  ultraviolet  radiation  bleaches  many 
fabrics.  Unless  a  sizeable  amount  of  heat  absorbing  mass  is 
present,  the  living  area  will  heat  up  quickly  during  daylight 
hours,  possibly  overheating.  Then,  during  the  night,  the 
building  will  lose  heat  through  the  glazing  by  convection 
and  radiation  from  the  building  interior  to  the  inner  surface 
of  the  glazing,  conduction  through  the  glazing,  and  convec¬ 
tion  and  radiation  from  the  exterior  surface. 

The  second  passive  heating  design  is  indirect  gain. 

Here,  heat  is  collected  and  transferred  first  to  the  thermal 
storage  mass.  Then,  the  heat  is  conducted  through  the  mass 
and  convected  and  radiated  to  the  interior.  The  most  common 
indirect  gain  system  is  the  Trombe  wall.  It  consists  of  a 
double  glazed  masonry  wall  that  has  the  glazed  side  exposed 
to  the  sun  and  the  other  side  facing  the  interior  of  the 
structure.  See  Fig  2.  Sunlight  warms  the  outside  surface 
of  the  wall  and  the  air  in  the  space  between  the  glazing 
and  the  wall.  Heated  air  from  the  air  space  can  be  circu¬ 
lated  through  vents  to  the  interior  to  obtain  immediate 
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energy  gains,  while  heat  conducted  through  the  wall  enters 

the  interior  some  hours  later. 

The  Trombe  wall  has  advantages  over  a  direct  gain 
system.  First,  glare  and  bleaching  are  not  a  problem.  Next, 
the  interior  is  not  subject  to  a  large  temperature  swing. 
Because  objects  inside  the  living  area  are  not  directly 
irradiated,  there  is  less  chance  of  overheating  during  the 
day.  Also,  because  the  interior  is  not  exposed  to  a  large 
expanse  of  glass,  heat  losses  are  held  down  at  night.  In 
fact,  a  Trombe  wall  will  continue  to  heat  room  air  long 
after  the  sun  has  gone  down.  For  a  twelve  inch  thick  ma¬ 
sonry  wall,  the  delay  from  the  time  the  sun  heats  the  out¬ 
side  of  the  wall  until  the  time  when  this  heat  soaks  through 
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to  the  Interior  can  be  eight  hours  or  more.  If  the  sun  sets 
at  4i 00  p.m. ,  this  Trombe  wall  will  continue  to  heat  the  in¬ 
terior  until  midnight. 

The  chief  disadvantage  of  the  Trombe  wall  is  its  heat 
loss.  The  most  efficient  way  of  heating  the  contents  of  the 
interior  during  the  day  is  by  directly  irradiating  them  with 
sunlight.  By  irradiating  the  outside  surface  of  a  wall, 
some  heat  soaks  into  the  wall,  but  loss  does  occur  through 
the  glazing  to  the  outside  in  the  same  fashion  as  through 
direct  gain  glazing.  This  loss  continues  at  night  as  well. 

If  the  storage  mass  could  be  insulated  from  the  out¬ 
side,  more  heat  would  find  its  way  to  the  interior,  and 
less  would  be  lost  to  the  outside.  One  way  to  accomplish 
this  is  to  detach  the  storage  mass  from  the  collector  and 
place  it  at  the  back  of  or  below  the  structure.  Heated  air 
from  the  collector  would  be  transferred  by  natural  convec¬ 
tion  or  by  fan  to  the  detached  storage  mass,  through  the 
mass,  and  back  to  the  a unspace.  See  Pig  3.  Because  the 
mass  is  well  insulated  from  the  outside,  more  heat  will  be 
available  for  the  interior.  In  addition,  the  south  side  of 
the  living  space  would  be  insulated  conventionally,  further 
reducing  heat  loss. 

When  designing  a  passive  solar  structure,  a  mix  of 
direct  and  indirect  features  might  be  incorporated.  The 
designer  would  try  to  maximize  the  average  temperature  gain 
over  outside  temperature,  keep  the  temperature  swing  within 
acceptable  limits,  and  minimise  the  amount  of  oostly  thermal 
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mass  that  needs  to  he  installed. 

Heat  Gain  In  the  Structure 

The  interior  of  a  structure  is  heated  by  the  sun  dir¬ 
ectly  and  indirectly*  by  internal  equipment  such  as  ovens 
and  water  heater? ,  by  people*  and  by  a  furnace.  Each  of 
these  mechanisms  for  heat  gain  will  be  discussed. 

Solar  Energy.  Just  outside  Earth's  atmosphere*  the 
energy  incident  on  a  flat  surface  directly  facing  the  sun 
would  be  13^2  M/m2  (432  B/f t2 ) .  As  sunlight  enters  the 
atmosphere*  the  rays  are  subjected  to  the  effects  of  the 
atmosphere.  Some  rays  are  reflected  back  by  clouds  into 
spaoe,  some  rays  art  scattered  or  absorbed  by  particles  or 
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air  molecules,  and  others  are  transmitted  to  the  ground 
where  they  are  absorbed  or  reflected.  The  result  of  these 
atmospheric  processes  is  that  an  object  on  the  ground  will 
receive  radiation  from  the  sun  in  three  different  forms. 

The  first  form  of  radiation  is  beam  radiation.  By 
definition,  this  is  radiation  where  the  light  rays  are 
very  nearly  parallel.  Beam  radiation  comes  directly  from 
the  sun  in  a  nearly  straight  line  (bent  only  a  slight 
amount  by  refraction)  to  the  object  in  question*  To  re¬ 
ceive  the  greatest  benefit  from  beam  radiation,  the  surface 
of  a  flat  collector  should  be  normal  to  a  line  drawn  between 
the  collector  and  the  sun. 

The  second  kind  of  solar  radiation  incident  upon  an 
object  is  diffuse  radiation  from  the  Bky.  Dust,  water  va¬ 
por,  pollution,  clouds,  and  the  air  itBelf  tends  to  scatter 
light.  The  result  is  that  light  reaches  an  object  on  the 
surface  from  all  directions  in  the  sky.  The  best  orienta¬ 
tion  for  a  flat  collector  to  receive  diffuse  sky  radiation 
is  face  up,  so  as  to  be  able  to  view  all  of  the  sky. 

The  last  type  of  solar  radiation  is  ground  radiation. 
Beam  and  diffuse  sky  radiation  reach  other  objects  and  the 
ground  just  as  it  reaches  the  collector.  Some  of  this  ra¬ 
diation  reflects  off  these  surfaces  toward  the  collector. 
Objects  such  as  glass,  water,  and  ice  can  reflect  specular¬ 
ly.  This  means  that  beam  radiation  is  reflected  as  beam 
radiation.  Objects  made  out  of  concrete,  grass,  and  earth 
reflect  diffusely,  which  means  that  incident  light  is 
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scattered  in  all  directions.  Ground  radiation  can  therefore 
he  either  beam  or  diffuse  in  nature .  However,  for  a  structure 
made  out  of  wood,  brick,  and  asphalt,  surrounded  by  grass, 
snow,  trees,  and  other  structures,  virtually  all  ground 
radiation  will  be  diffuse. 

The  amount  of  beam,  diffuse,  and  ground  radiation  that 
a  collector  receives  is  dependent  on  many  factors,  some 
of  which  were  mentioned  earlier.  Collector  orientation  is 
important.  A  collector  facing  straight  up  will  receive  no 
ground  radiation.  A  collector  looking  straight  north  will 
receive  little  beam  radiation,  for,  in  the  Northern  Hemis¬ 
phere,  the  sun  is  rarely  found  north  of  an  east-west  meri¬ 
dian.  Another  factor  that  affects  insolation  is  the  sky 
condition.  If  the  sun  is  obscured  by  clouds,  no  beam 
radiation  will  reach  the  collector.  Also,  the  ratio  of 
diffuse  to  beam  radiation  will  be  greater  in  New  Orleans 
on  a  humid  day  than  in  Albuquerque  on  a  dry  day.  Other 
factors  that  affect  the  kind  and  amount  of  radiation  reach¬ 
ing  a  collector  include  time  of  day,  building  overhangs, 
adjacent  trees  or  structures  blocking  morning  or  evening 
sunlight,  the  location  of  the  collector,  and  the  time  of 
year. 

Over  the  years,  solar  radiation  data  collection  sta¬ 
tions  have  been  set  up  around  the  world.  Information  on 
hourly,  daily,  monthly,  and  annual  total  radiation!  clear¬ 
ness  indiciesi  and  temperatures  are  available  in  tabular 
form.  Perhaps  the  most  readily  available  data  is  monthly- 
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average  daily  total  radiation  on  a  horizontal  surface,  monthly- 
average  clearness  index,  and  monthly-average  temperature. 

To  be  useful  for  simulations,  solar  insolation  values 

« 

for  any  time  during  the  day  must  be  available.  Duffie  and 
Beckman  (Ref  4)  devised  an  algorithm  for  converting  monthly- 
average  daily  total  radiation  into  beam  and  diffuse  compon¬ 
ents  of  hourly  radiation.  Before  outlining  the  technique, 
a  few  terms  need  to  be  defined. 

The  hour  angle  CJ  is  the  angular  displacement  of  the  sun 
east  or  west  of  the  local  meridian  due  to  the  rotation  of 
the  earth  on  its  axis,  i.e.,  15  degrees  per  hour.  The  hour 
angle  corresponding  to  solar  noon  is  zero.  Morning  hour 
angles  are  negative.  The  hour  angle  is  given  in  degrees  by 
CJ  =  15  (TIME  OF  DAY)  -  180 
where,  for  example,  1330  hours  would  be  13*5. 

The  declination  S  is  the  angular  position  of  the  sun 
at  solar  noon  with  respect  to  the  plane  of  the  equator. 

Because  the  earth's  axis  is  tilted  23 .45  degrees  to  the 
normal  to  the  plane  of  the  ecliptic,  the  declination  will, 
depending  on  the  day  of  the  year,  be  at,  or  within  limits 
of  +  23.45  degrees. 

S  =  23.45  sin(360(284  +  (DAY  OF  YEAR))  /  365) 

The  solar  collector  can  be  located  at  latitude  <f)  ,  and 
be  oriented  as  to  azimuth  ,  and  tilt  /?.  If  a  collector 
is  facing  south,  it  has  an  azimuth  of  zero.  If  it  is  facing 


13 


east  of  south,  its  azimuth  is  negative.  If  a  collector  is 
horizontal,  its  tilt  is  zero.  If  it  is  vertical,  the  tilt 
is  90  degrees. 

The  sunset  hour  angle  is  the  hour  angle  at  which  the 
sun  sets. 

CJQ  =  arceos(-tan  <f>  tan 8  ) 

Given  the  monthly-average  daily  total  radiation  H  and 
the  monthly-average  clearness  index  K^,  the  monthly-average 
daily  diffuse  radiation  Bd  is  calculated  with  the  following 
relation! 

Hd  =  H  (0.775  +  0.00653  (CJB-  90)  -  (0.505  + 

0.00455  ((*)B  -  90))  cos(115  Kt  -  103))  (1) 

The  hourly  totul  radiation  I  is 

I  =  H  ((7 T/  24)  (a  +  b  cosCJ)  (cos  CJ  -  cosCJ8)  / 

(sin£Js  -  (2  7TCJb  /  360)  cos  U)B))  (2) 

where 

a  =  0.409  +  0.5016  sin( CO  -  60) 

D 

and 

b  =  0.6609  -  0.4767  sin(CJB  -  60) 

The  hourly  diffuse  radiation  Id  is 

Id  *  ((7T/  24)  (cos  U)  -  cos  CJ8)  / 

(sin  CJ8  -  (2  7T  /  360)  cos  CJ8)) 
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Then,  hourly  beam  radiation  iB  the  hourly  total  less  the 
hourly  diffuse  radiation. 

h  -  1  -  *d  (*> 

The  values  of  hourly  beam  and  diffuse  radiation  calcu¬ 
lated  with  the  above  relations  are  for  a  collector  that  is 
facing  straight  up.  The  collector  for  a  residential  size 
structure  passively  heated  by  the  sun  will,  in  the  vast 
majority  of  the  cases,  be  a  glass  window  hung  vertically, 
facing  generally  south.  The  values  for  hourly  beam  and  dif¬ 
fuse  radiation  must  be  altered  for  a  collector  that  has  this 
orientation. 

Hourly  beam  radiation  1^  must  be  multiplied  by  a  geo¬ 
metric  factor  R^  which  is  the  ratio  of  beam  radiation  on 
a  tilted  surface  to  that  on  a  horizontal  surface. 

=  cos  6  /  cos  ez  (5) 

where  the  angle  of  incidence  8  is  the  angle  between  the 
direction  of  beam  radiation  and  the  normal  to  the  collector 
surface.  The  zenith  angle  6Z  is  the  angle  between  the  direc¬ 
tion  of  beam  radiation  and  the  normal  to  a  horizontal  surface. 
The  angle  of  incidence  and  the  zenith  angle  can  be  calculated 
from  the  following  relations t 

cos  6  =  sin  8  si n(f)  cos  ft  -  sin  8  coB<f>  sin  ft  cos  Y  + 
cos  8  coeCpaoB  ft  oob  U)  +  cos  5  sin  <f>  sin  ft 
cos  y  cob  U)  +  cos  8  sin  ft  sin  y  sin  OJ 
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cos  6Z  -  sin  8  sin  (f)  +  cos  8  cos  <fi  cos  CJ 

For  -a  collector  oriented  vertically,  the  tilt  will  be 
90  degrees.  The  equation  for  the  angle  of  incidence  is 
simplified  to 

cos  0  =  -sin  8  cos  <f>  cos  y  +  cos  8  sin  <f)  cos  y  cos  CJ 
+  cos  8  sin  *Y  s^n  CJ  (7) 

Diffuse  radiation  has  no  preferred  direction.  ItB  in¬ 
tensity  is  the  same  no  matter  which  way  the  collector  is 
oriented.  To  maximize  the  collection  of  diffuse  radiation, 
the  collector  should  be  exposed  to  as  much  of  the  sky  as 
possible.  For  a  window  in  the  wall  of  a  house,  only  half 
the  sky  can  be  sen.  Overhangs,  trees,  and  adjacent  build¬ 
ings  will  cut  out  additional  portions  of  the  sky.  For  this 
study,  it  will  be  assumed  that  one  half  of  the  sky  can  be 
seen  by  the  collector  surface.  Therefore,  half  of  the 
hourly  diffuse  radiation  1/2  Id  is  incident  on  the 
collector. 

Additional  radiation  is  incident  on  the  vertically 
oriented  collector.  The  amount  of  hourly  ground  radiation 
incident  on  a  flat  surface  facing  straight  down  is  the 
product  of  the  hourly  total  radiation  I  and  the  reflec¬ 
tivity  of  the  ground  /^.  Again,  as  the  collector  is  verti¬ 
cal,  and  the  ground  radiation  is  diffuse,  only  half  of  the 
ground  radiation  available  will  be  incident  on  the  collector 
surface.  This  relation  is  given  by  1/2  . 
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Given  the  solar  radiation  incident  on  the  collector 
surface*  the  rate  of  heat  gain  for  the  structure  can  be 
calculated.  For  a  direct  gain  system,  the  rate  of  solar 
gain  Qpg  across  a  collector  of  area  can  be  written  as 

Ono  -  Ajo  O'  <ib  Rb  Tb  +  1/2  i„  Td  +  1/2  Pg  I  V  (8) 

where  &  is  the  effective  absorptance  of  the  window-room 
combination,  and  7^  and  7^  are  the  beam  and  diffuse  trans- 
mittances  of  the  window  glazing,  respectively.  The  com¬ 
bination  absorptance  OC can  be  approximated  by 

...  .  ^d  Pr  ADG  /  AR 

where  PR  is  the  reflectance  of  the  interior  walls,  floor 
and  ceiling,  and  AR  is  the  area  of  the  interior  walls, 
floor  and  ceiling. 

The  heat  collected  each  hour  at  the  outer  surface  of 
a  Trombe  wall  of  area  A^  can  be  expressed  as 

QTO  *  ATW  *xb  Rb  ^b  +  1//Z  *d  ^d  +  /%  1  *^d(9) 

where  70^  and  70^  are  the  effective  beam  and  diffuse 
transmittance  through  the  glazing  and  absorptance  of  the 
masonry  wall,  respectively.  This  heat  warms  the  air  in  the 
space  between  the  glazing  and  the  wall,  and  conducts  into 
the  wall  itself  where  it  is  stored.  While  the  air  from 
the  Bpac e  may  be  circulated  through  vents  into  the  interior 
to  affect  immediate  interior  temperature  gains,  this  option 
will  be  omitted  in  this  analysis  for  two  reasons.  First, 


heat  that  is  immediately  circulated  cannot  he  stored.  This 
reduces  the  effectiveness  of  the  delayed  action  heating 
that  makes  the  Trombe  wall  so  popular.  Second,  a  combina¬ 
tion  of  direct  gain  heating  for  early  morning  gains  and 
indirect  gain  for  later  heating  makes  the  option  unnecessary. 

For  a  detached  storage  mass,  one  must  calculate  the 
heating  of  the  air  in  the  sunspace  and  the  transport  of 
heat  to  the  storage  mass.  Heat  gain  by  any  masonry  in  the 
sunspace  can  be  calculated  using  Eq  9.  substituting  the 
surface  area  of  the  mass  in  the  sunspace  A«,s  exposed  to  the 
sky  for  the  surface  area  of  the  Trombe  wall  A^.  This  mass 
will  heat  the  air  in  the  sunspace  by  convection. 

=  ASS  ^  "  TSS*  (10> 

where  QgS  is  the  rate  of  heat  transfer  from  the  mass  in  the 
sunspace  to  the  air,  hgS  is  the  convection  coefficient,  T^ 
is  the  surface  temperature  of  the  mass  in  the  Bunspace,  and 
Tss  is  the  temperature  of  the  air  in  the  sunspace.  Heated 
air  from  the  collector  is  transported  to  the  storage  mass 
by  natural  or  forced  convection,  then  returned  to  the  sun¬ 
space  for  further  heating.  The  rate  of  energy  transport  to 
the  detached  storage  mass  is  A  Cp  TSS’  where  A  is  the  mass 
flow  rate  and  Cp  is  the  specific  heat  of  air.  The  mass 
flow  rate  is,  if  by  natural  convection,  dependent  on  the 
geometry  of  the  system  and  the  temperature  gradients.  Mass 
flow  rates  for  foroed  convection  depend  on  the  size  of  the 
fans  and  resistance  to  flow  in  the  ductwork.  If  the  flow  is 
forced,  the  sunspace  air  oould  be  pushed  through  an  efficient 
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heat  exchanger  until  it  reaches  thermal  equilibrium  with 
the  center  of  the  storage  mass.  If  this  is  the  case*  the 
rate  at  which  heat  is  transferred  to  the  mass  is 
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where  is  the  heat  flow  rate  and  is  the  temperature 

of  the  center  of  the  storage  mass. 

Heat  Gain  from  Won-Solar  Sources.  The  heat  generated 
by  appliances  and  people  within  a  structure  is  not  insig¬ 
nificant.  Schools,  shopping  centers,  and  office  buildings 
often  times  have  overheating  problems  with  freezing  tempera¬ 
tures  outside.  As  an  example,  an  average  human  body  lying 
down  will  generate  300  B/hr.  Tiis  heat  reduces  the  ad¬ 
ditional  energy  that  must  be  supplied  by  the  furnace  or 
the  solar  heating  system. 

A  residential  structure,  with  few  people  and  little 
equipment,  has  a  small  amount  of  internal  heat  generated. 
According  to  the  ASHRAE  Handbook  (Ref  5)»  the  internal  sources 
in  a  residential  structure  are  good  for  a  7  degree  Fahrenheit 
increase  in  the  interior  temperature.  Accepted  practice  in 
the  industry  is  to  disregard  all  internally  generated  heat 
and  design  the  heating  system  to  maintain  an  interior 
temperature  of  65  F. 

The  furnace  in  a  residential  structure  is.  most  often, 
an  on-off  heat  generator  that  heats  the  interior  air  at  a 
constant  rate  when  on.  The  furnace  is  controlled  by  a 
thermostat,  which  turns  the  furnace  on  when  the  interior 
temperature  reaches  a  lower  limit  and  turns  it  off  when  an 


upper  temperature  limit  is  reached.  In  this  study,  the 
lower  limit  will  be  set  at  65  F,  the  upper  limit  at  69  F. 

Heat  Loss  From  the  Structure 

During  the  winter  months*  valuable  heat  from  the 
interior  of  a  structure  finds  its  way  to  the  outside  and 
is  lost.  Heat  is  lost  by  conduction  through  the  walls* 
ceiling,  and  roof,  then  convection  and  radiation  to  the  out¬ 
side  air.  Heat  is  also  lost  by  conduction  through  the  ground. 
Lastly*  heat  is  lost  by  infiltration  of  outside  air  in  and 
inside  air  out.  Each  of  these  losses  will  be  examined. 

Loss  through  Walls.  Windows,  and  Ceiling.  This  loss 
is  calculated  using  the  relation  Q  =At  /  R  .  where  Q  is 
the  rate  of  heat  Ioss.At  is  the  difference  between  interior 
and  outside  temperature*  and  R  is  the  total  thermal  resistance. 
The  resistance  is  calculated  by  summing  the  thermal  resist¬ 
ances  of  the  components  of  the  wall,  ceiling,  or  window. 

Heat  loss  through  the  foundation  is  a  different  matter. 

The  ground  is  a  good  insulator.  Little  heat  is  lost  straight 
down.  Most  of  the  floor  loss  is  around  the  perimeter*  where 
the  distance  to  the  outside  air  is  at  its  shortest.  Most 
new  houses  built  in  cold  climates  have  perimeter  insulation. 

The  resistance  to  heat  loss  in  the  perimeter  is  a  function 
of  its  length*  the  thickness  of  the  insulation,  and  the  depth 
below  ground  level  to  which  the  insulation  extends.  For 
english  units*  perimeter  resistance  is  given  by 
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Table  I  Edge  Loss  Factor  F  for  a  Masonry  Slab  Floor 
Insulation  R-Value  Depth  of  Insulation 


1  ft 

1.5  ft 

2  ft 

10.0 

0.20 

0.18 

0.17 

6.7 

0.32 

0.28 

0.27 

5.0 

0.43 

0.39 

0.37 

4.0 

0.57 

0.50 

0.48 

3.3 

0.70 

0.62 

0.59 

2.5 

1.00 

0.88 

0.83 

From t  Mazria  (Ref  6) 

where  P  is  the  length  of  the  perimeter  and  F  is  an  edge 
loss  factor  from  Table  I. 

Perhaps  the  largest  single  factor  contributing  to 
heat  loss  is  infiltration.  Outside  air  seeps  in,  and  in¬ 
side  air  out  through  cracks  around  doors  and  windows, 
through  stove  and  dryer  vents  and  furnace  flues,  and 
through  the  building  materials  themselves.  In  addition, 
every  time  a  door  or  window  is  opened,  additional  air  is 
exchanged.  For  a  building  with  no  weather  stripping  or 
storm  windows,  the  infiltration  rate  can  be  as  high  as 
four  complete  air  changes  per  hour.  By  weatherstripping 
and  sealing  all  cracks,  this  rate  can  be  brought  down  to 
about  3/4  change  per  hour.  Even  at  this  rate,  though, 
infiltration  can  account  for  40  percent  of  the  total 
building  heating  load.  Further  reductions  would  necessitate 
installing  vapor  seals,  either  paint  or  plastic  sheeting  on 
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the  walls  and  ceiling.  However,  fresh  air  is  important  to 
the  well-being  of  the  occupants.  Below  one  half  an  air 
change  each  hour,  carbon  dioxide,  water  vapor,  and  cooking 
odors  build  up.  The  resistance  to  heat  loss  is  given  by 
the  relation 


where  is  the  specific  heat  of  air,  /3  is  the  density  of 
air,  V  is  the  volume  of  the  structure,  and  n  is  the  exchange 
rate. 

Heat  Capacitance  in  the  Structure.  The  materials  that 
go  into  a  building,  wood,  bricks,  concrete,  all  have  a 
capacity  to  store  heat.  This  thermal  mass  heats  up  and 
cools  off  more  slowly  than  air  surrounding  the  mass.  A 
large  amount  of  mass  will  moderate  the  temperature  swing 
inside  the  structure.  Calculating  the  heat  capacity  in¬ 
volves  finding  the  product  of  the  mass  of  each  building 
component  and  its  specific  heat.  For  example,  a  20  ft  by 
8  ft  by  1  ft  thick  brick  wall  weighing  19200  pounds  and 
having  a  specific  heat  of  0.2  B/lbm-F  will  have  a  heat 
capacity  of  3840  B/F.  This  means  that  for  the  wall  to 
change  one  degree  in  temperature,  it  must  gain  or  lose 
3840  B  from  or  to  its  surroundings. 

Economics 

Despite  the  advantages  that  solar  heating  has,  i.e., 
self-sufficiency  and  thermal  comfort,  very  few  solar  homes 


would  be  built  if  it  were  not  economical  to  do  so.  The 
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major  disadvantage  of  building  a  solar  structure  ie  its  high 
initial  cost.  Passive  solar  structures  require  a  ggreat  deal 
of  glazing  with  its  necessary  framing*  and  heavy  walls  and 
floors.  While  frames  and  walls  can  be  incorporated  into 
the  support  structure  of  the  building,  the  construction 
cost  for  a  solar  dwelling  is  higher  than  that  for  a  conven¬ 
tional  one.  On  the  other  hand*  a  well  designed  passive 
solar  home  will  reduce  monthly  heating  bills,  and  initial 
cost  can  be  partially  offset  by  tax  incentives.  Federal 
and  state  governments  allow  the  builder  who  uses  solar 
energy  to  deduct  some  construction  costs  from  taxable  in¬ 
come. 

Life  Cycle  Savings  Analysis  is  useful  in  determining 
the  advisability  of  building  a  solar  dwelling.  First,  the 
present  worth  of  additional  costs  for  construction,  in¬ 
surance,  property  taxes  and  mortgage  payments  is  computed, 
taking  into  account  inflation  and  interest  rates,  over  a 
projected  life  span,  typically  20  or  30  years.  Then  the 
same  calculations  are  made  for  utility  savings,  tax  credits, 
and  tax  deductions  for  additional  mortgage  interest  payments. 
The  Solar  Savings  is  obtained  by  subtracting  the  costs  from 
the  savings.  A  negative  result  indicates  that  the  solar 
house  will  be  more  expensive  than  a  conventional  one.  A 
positive  result  means  that,  in  the  long  run,  the  solar 
house  will  save  the  owner  money. 


Thermal  Comfort 


-  What  a  person  perceives  as  the  temperature  in  a 
dwelling  is  actually  a  function  of  several  factors.  The 
temperature  of  the  air  immediately  surrounding  a  body  and 
the  velocity  of  the  air  past  the  body  affect  this  perception. 
Cooler  air,  and  moving  air  both  increase  the  rate  of  heat 
transfer  from  a  warm  body  to  its  surroundings.  In  addition, 
the  humidity  of  the  air  is  a  factor.  Moisture  in  the  air 
decreases  the  rate  of  evaporation  of  sweat,  decreasing  the 
heat  loss  by  transpiration  cooling.  Also,  the  activity 
level  and  the  amount  of  clothing  a  person  wears  affects  how 
warm  that  person  feels.  The  last  and  perhaps  most  neglected 
factor  is  the  temperature  of  the  walls,  floor  and  ceiling. 
Thermal  energy  radiated  from  these  surfaces  Btrikes  and 
warms  a  person.  The  absence  of  radiation  from  a  warm  wall 
can  easily  be  detected.  Standing  in  a  warm  room  next  to  a 
picture  window  during  the  winter,  a  person  will  feel  colder 
on  the  side  that  is  facing  the  window,  for  less  radiant 
energy  is  incident  on  a  body  from  the  cold  glass  than  from 
a  warm  wall. 

Wray  (Ref  ?)  outlines  a  method  for  determining  thermal 
comfort  which  takes  into  account  all  of  these  factors.  In 
a  dwelling,  the  velocity  of  air  is  not  sufficient  to  dis¬ 
rupt  natural  convection.  Therefore,  motion  of  air  can  be 
neglected.  In  addition,  for  reasons  of  comfort,  the  rela¬ 
tive  humidity  must  be  at  moderate  levels.  According 

to  Wray,  comfortable  equivalent  uniform  temperatures  for 
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relative  humidities  of  zero  and  fifty  percent  are  shown  in 
Table  II.  The  equivalent  uniform  temperature  Tgu  is  the 
uniform  temperature  of  an  imaginary  enclosure  in  which  a 
person  will  experience  the  same  degree  of  thermal  comfort 
as  in  the  actual  non-uniform  environment.  For  a  relative 
humidity  of  zero  percent, 

Teu  =  °.52  T^.  +  0.48  Ta  (14) 

where  Tmr  is  the  mean  radiant  temperature  of  the  walls, 
floor  and  ceiling,  and  Tj.  is  the  interior  air  temperature. 

For  a  relative  humidity  of  fifty  percent, 

T„u  -  Tmr  +  "•»  Ti  (15> 

The  mean  radiant  temperature  Tmr  of  the  walls,  ceiling 
and  floor  is  the  hypothetical  uniform  temperature  of  the 
interior  surfaces  that  would  result  in  the  same  heat  trans¬ 
fer  due  to  radiation  to  an  object  in  the  interior  of  the  room. 


Table  III  Comfort  Range  At 
0%  Rel.  Humidity  50%  Rel.  Humidity 

Activity  Level  Activity  Level 


Clothing 

Office 

Work 

Domestic 

Work 

Office 

Work 

Domestic 

Work 

Light  slacks, 
short  sleeve 
shirt 

5.5  F 

3  C 

10  F 

5*5  C 

4.5  F 

2.5  c 

9  F 

5  C 

Typical 

business 

8  F 

4.5  C 

14.5  F 

8  C 

5.5  F 

3  C 

12.5  F 
7  C 

suit 

Extracted  fromt  Wray 


The  temperature  comfort  range  centers  on  the  equiva¬ 
lent  uniform  temperature  Tgu.  According  to  Wray,  a  human 
subject  will  experience  no  more  than  slight  discomfort  when 
the  equivalent  uniform  temperature  is  in  the  temperature 
range 

T  =  T  ,  +  At  (16) 

r  eu  - 

where  T  is  the  comfort  range  from  Table  III. 

Optimization 

Comfort  and  cost  are  the  overriding  factors  to  be 
considered  in  selecting  a  design  for  a  dwelling.  No  house 
that  is  uncomfortable  is  worth  living  in.  However,  given 
that  the  building  is  designed  so  an  adequate  average  temp¬ 
era  ^ ture  is  maintained,  and  that  the  temperature  swing  is 
held  within  tolerances  described  earlier,  the  cheapest 
design  will  be  the  best. 


Ill  The  Computer  Model 

The  approach  used  to  analyse  the  heat  flow  in  a  passive 
solar  structure  was  to  construct  an  electrical  analog  circuit. 
The  thermal  resistance  of  the  building  components  were  model¬ 
ed  with  electrical  resistors.  Energy  storage  capability  was 
simulated  with  capacitors.  Heat,  then,  was  the  flow  between 
adjacent  structural  components  with  different  temperatures 
just  as  current  flows  between  adjacent  nodes  with  different 
voltages . 

The  Circuit  Diagram 

The  circuit  diagram  for  the  structure  is  shown  in 
Fig  4.  By  analysing  the  circuit  with  different  values  for 
the  resistors  and  capacitors,  and  by  specifying  different 
values  for  solar  insolation  and  outside  temperature,  simu¬ 
lations  for  a  wide  variety  of  structure  configurations  and 
sizes  can  be  run  for  any  location  across  the  country  and 
any  month  of  the  year.  An  explanation  of  each  component  of 
the  circuit  and  its  real  world  analog  will  follow. 

The  Ambient  Temperature  TQ.  Data  available  in  solar 
insolation  tables  include,  for  each  location,  average  out¬ 
door  temperatures,  including  average  highs  and  lows,  for 
each  month  of  the  year.  For  purposes  of  the  simulation, 
ambient  temperature  was  approximated  by 

Ta  *  Tave  “  Tswing  ^ob((TIUE  OF  DAY  -  2)  7 T/  12)  (17) 

where  Tftye  is  the  mean  monthly  temperature  and  Tgw^n^  is 
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Fig  5.  Schematic  and  Circuit  Diagram  for  Trombe  Wall 


half  the  difference  between  the  average  high  and  the 
average  low  for  that  month. 

Nodes  1  through  Trombe  Wall.  Refer  to  Fig  5. 

Given  that  the  surface  area  is  large  compared  to  the 
thickness,  one  dimensional  heat  flow  is  a  good  assump¬ 
tion.  According  to  Duffie  and  Beckman,  heat  flow  in  a 
thick  wall  can  be  adequately  modeled  with  five  nodes.  Fol¬ 
lowing  this  suggestion,  one  node  was  placed  at  the  outside 
surface,  one  at  the  inside  surface,  and  three  equally  spaced 
within.  Total  heat  capacitance  of  the  wall  was  C p  P  V  , 
where  Cp  is  the  specific  heat  of  the  masonry,  p  is  its 
density,  and  V  is  its  volume.  Capacitors  C2,  C3,  and  C4 
were  each  assigned  values  of  one  fourth  the  total  wall 
capacitance,  while  capacitors  Cl  and  C5  were  each  assigned 
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one  eighth  the  total  value.  Total  internal  resistance  was 
given  by  k  V  /  A,^  ,  where  k  is  the  resistance  to  conduc¬ 

tion  in  the  masonry,  and  is  the  cross  sectional  area  of 
the  wall.  Resistors  R3  through  R6  were  each  assigned  a  value 
of  one  fourth  the  total  internal  resistance. 

Resistor  R2  models  the  combined  resistance  of  the 
natural  convection  air  film  between  the  sunspace  and  the 
glazing,  resistance  to  conduction  in  the  glazing,  and  insula¬ 
tion  provided  by  the  air  gaps.  Resistor  R?  represents  the 
natural  convection  air  film  between  the  wall  and  the  inter¬ 
ior  air  space. 

Solar  input  for  the  Trombe  wall  is  added  at  node 
1.  The  value  of  this  input  was  calculated  using  the  algo¬ 
rithm  detailed  in  Section  II, 

Nodes  6  to  10 t  Detached  Storage  Mass.  As  shown  in  Fig 
6,  the  model  for  the  detached  storage  mass  is  very  similar 
to  the  model  for  the  Trombe  walli  five  nodes  spanning  the 
thinnest  dimension  (the  direction  of  heat  flow).  Values  for 
capacitors  C6  through  CIO  and  resistors  Rll  through  R14  are 
calculated  the  same  way  as  for  capacitors  Cl  through  C5  and 
resistors  R3  through  R6,  respectively.  For  this  model, 
heated  air  from  the  sunspace  passes  down  channels  through 
the  center  of  the  storage  mass.  A  good  configuration  would 
be  blowing  air  through  a  cinder  block  wall  with  the  cavi¬ 
ties  aligned.  For  this  reason,  resistor  R9  connects  the 
sunspace  with  node  8.  A  method  for  calculating  the  value 
of  resistor  R9  is  shown  in  Section  II.  Resistor  RIO  may  be 
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the  perimeter  resistance  if  the  detached  storage  mass  is 
in  the  floor  of  the  structure)  it  is  the  resistance  of  an 
outside  wall  if  the  storage  mass  is  at  the  back  (north  side) 
of  the  structure. 

Nodes  11  and  12 >  Outside  Walls  and  Roof/Celling. 

The  models  for  a  wall  and  a  roof/ceiling  with  insulation 
are  shown  in  Fig  7.  Note  that  mass  outside  the  insulation 
envelope  is  not  assigned  a  capacitance.  Since  structural 
components  outside  the  insulation  are  not  thermally  pro¬ 
tected  from  outside  temperatures,  these  temperatures  were 
considered  to  be  close  to  ambient  temperature.  Materials 
inside  the  installed  insulation  were  considered  to  be 
thermal  mass  affecting  the  interior  temperature  swing,  and 
were  considered  when  calculating  the  values  of  capacitors 
Cll  and  C12.  Values  for  resistors  R16  and  R18  were  obtained 
by  summing  the  thermal  resistances  for  the  building  compon¬ 
ents  involved.  Resistors  R17  and  R19  modeled  resistance 
due  to  convection  air  films  on  the  walls  and  ceiling,  re¬ 
spectively.  The  area  for  determining  the  wall  capacitance 
and  resistances  is  the  vertical  area  (product  of  perimeter 
length  and  wall  height)  less  all  window  area  and  wall  area 
accounted  for  by  the  Trombe  wall  and  or  the  detached  stor¬ 
age  mass.  The  area  for  determining  the  roof/ceiling  capaci¬ 
tance  and  resistances  is  the  interior  floor  area. 

Nodes  13.  14.  and  17 i  Floor.  This  model  allows  for 
a  masonry  or  concrete  slab  floor,  perhaps  with  tile  or 
linoleum  covering.  The  floor,  being  thinner  than  a  Trombe 
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.  '  ,  wall  or  detached  storage  mass,  can  be  modeled  with  fewer 

than  five  nodes.  Three  nodes  were  selected,  as  shown  in 
Fig  8.  Values  for  capacitors  C13,  Cl4,  and  C15  and  resis¬ 
tors  R20  and  R21  are  determined  in  the  same  manner  as  for 
the  Trombe  wall.  Capacitor  C13  has  a  value  of  one  half  the 
total  capacitance,  while  capacitors  C14  and  C15  each  have 
a  value  of  one  fourth.  Resistors  R20  and  R21  are  each 
assigned  values  of  one  half  the  internal  resistance  of  the 
floor.  Resistor  R24  is  an  air  film  resistance  between  the 
floor  and  the  interior,  while  resistor  R25  is  the  perimeter 
resistance  described  in  Section  II. 

The  direct  gain  insolation  that  enters  the  structure 
is  absorbed  and  or  reflected  by  surfaces  in  the  room.  Walls, 

( 
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ceiling,  furniture,  and  floor  all  absorb  some  portion  of 
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Fig.  9  Circuit  Diagram  for  Heat  Flow  in 
S unspace  and  Interior 


direct  gain  radiation.  However,  for  most  structures,  in 
most  locations,  for  the  longest  portion  of  the  insolation 
period,  most  of  the  direct  radiation  is  incident  on  the 
floor.  For  this  reason,  the  direct  heat  gain  Q^g,  which 
is  calculated  using  the  method  described  in  Section  II, 
is  input  at  node  14. 

Nodes  IS  and  16 «  Interior  Air  and  the  Sunsoace  Air. 
These  do  not  have  capacitors  associated  with  them.  See  Fig 
9.  While  air  does  have  a  heat  capacitance,  it  is  small 
compared  to  that  of  the  thermal  mass.  Hence,  interior  air 
will  change  temperature  much  more  quickly  in  response  to  a 
given  heat  input  than  will  building  materials  in  the  struc¬ 
ture.  Resistors  Rl,  R8  and  R22  model  the  thermal  resistance 


in  the  glazing  between  outside  and  sunspace,  between  sunspace 
and  interior  through  the  direct  gain  glazing,  and  between  the 
interior  and  outside  through  non-south  windows,  respectively. 
Curtains,  night  insulation,  and  shades  can  affect  the  value 
of  these  resistors.  Resistor  R23  models  the  resistance  to 
infiltration  described  in  Section  II. 

The  interior  air  is  heated  by  convection  from  warm 
surfaces  such  as  the  interior  of  the  Trombe  wall  and  the 
floor.  Heat  gain  from  appliances  and  from  the  furnace 
QAUx  also  directly  heat  the  air.  While  heat  generated  by 
appliances  and  people  is  neglected  in  this  study,  the  pro¬ 
vision  for  specifying  internal  loads  is  included. 

Design  Choices 

Conventional  Structures  and  Direct  Gain  Systems.  Refer¬ 
ring  to  Fig  4,  conventional  structures  and  direct  gain  sys¬ 
tems  can  be  modeled  by  proper  selection  of  resistor  values. 

By  setting  resistor  R1  to  a  very  low  value,  say  10-1^,  resist¬ 
ors  R2.  R7.  R 9*  RIO,  and  R15  to  very  high  values,  say  10*^, 
and  setting  the  area  of  the  Trombe  wall  to  zero,  the  sunspace, 
Trombe  wall,  and  detached  storage  mass  are  effectively  elimin¬ 
ated  from  the  circui*.  Direct  gain  solar  heating  through  un¬ 
shuttered,  south  facing  glazing  can  then  be  modeled  by  making 
the  south  facing  glazing  area  small  for  a  conventional  house 
or  by  making  the  area  large  for  a  direct  gain  system. 

Indirect  Gain  Systems.  Again  referring  to  Fig  4,  in¬ 
direct  gain  systems  can  be  modeled  by  setting  resistors  R8, 


R9*  RIO,  and  R15  to  very  high  values,  and  setting  the  area 
of  the  south  facing  glazing  to  zero,  thereby  eliminating  the 
direct  gain  and  detached  storage  from  the  circuit.  The  re¬ 
maining  circuit  then  models  heat  flow  in  a  structure  with 
a  solar  heating  system  that  is  exclusively  Trombe  wall. 

To  simulate  a  structure  with  a  sunspace  and  detached 
storage  mass,  and  no  windows  between  the  sunspace  and  in¬ 
terior,  still  other  choices  are  required.  The  sunspace 
requires  mass  to  absorb  incident  radiation  and  heat  the 
contained  air.  Mass  in  the  form  of  a  vertical,  south 
facing  wall  can  be  modeled.  This  wall,  if  thin,  can  be 
modeled  with  node  1.  A  thicker  wall  in  the  sunspace  may 
require  use  of  nodes  2,  3  and  perhaps  4  to  yield  an  ade¬ 
quate  model.  Those  nodes  in  the  Trombe  wall  that  are  not 
used  to  model  the  sunspace  mass  must  be  isolated.  A  con¬ 
ventional  wall  used  to  separate  the  sunspace  from  the  in¬ 
terior  can  be  modeled  by  assigning  resistor  R8  a  value 
appropriate  to  simulate  resistance  to  heat  flow  in  such  a 
wall.  Lastly,  if  the  detached  storage  mass  is  also  to  be 
the  floor  of  the  structure,  resistors  R24  and  R25  must  be 
assigned  large  values  to  isolate  the  three-noded  model. 

Combination  of  Features.  With  imagination,  a  structure 
with  any  combination  of  the  above  passive  designs  may  be 
easily  modeled.  The  engineer  can  alter  his  or  her  design, 
or  alter  the  proportion  of  collector  dedicated  to  each 
type  of  system  to  find  the  most  effective  design  for  the 
location,  climate,  and  intended  use  of  structure. 

Wight  Configuration.  During  the  time  the  sun  is  not 


shining  on  the  collector,  heat  is  lost  through  it.  Because 
a  large  expanse  of  glass  has  very  low  resistance  to  heat 
flow,  many  builders  make  provisions  for  insulated  shutters, 
curtains,  or  shades  to  be  used  at  night.  For  this  reason, 
resistor  Rl  must  have  a  daytime  and  a  nighttime  value. 

In  addition,  air  flow  to  the  detached  storage  mass  may 
be  cut  off  at  night.  Therefore,  resistor  R9  also  must  have 
a  daytime  and  a  nighttime  value. 

The  Governing  Nodal  Equations 

Heat  flow  at  each  node  in  the  circuit  diagram  can  be 
described  by  a  governing  equation.  Consider  a  general  node 
x  attached  to  capacitor  x  and  connected  to  an  arbitrary  number 
n  of  nodes  as  shown  in  Fig  10.  A  general  heat  flow  Qx  is  also 
input  directly  to  the  node.  The  heat  flowing  to  node  x 
through  each  resistor  Rn  is  proportional  to  the  tempera tui*e 
difference  across  the  resistor  and  inversely  proportional 
to  the  magnitude  of  the  resistance.  This  is  given  by 


Tn  -  Tx 
Rn 


where  Tn  and  Tx  are  the  temperatures  at  nodes  n  and  x, 
respectively.  The  total  heat  flux  to  node  x  through  the 
resistors  and  as  a  result  of  the  direct  heat  input  Qx  is 

+  3^  +  +  ...  +  Tn^ 

The  rate  of  rise  in  temperature  at  node  x  is  proportional 
to  the  total  heat  flux  and  inversely  proportional  to  the 
value  of  capacitor  Cx. 

dTx  _  TOTAL  HEAT  FLUX 
W  *  Cx - 

where  T  is  time.  Hence,  the  governing  differential  equation 
for  node  x  is 

=  ^  +  +  T^c  +  ...  + 

If  the  value  of  capacitor  Cx  is  zero,  i.e.,  if  node  x  nas 
no  capacitor,  then  the  heat  flux  to  node  x  must  sum  to  zero. 
The  resulting  algebraic  equation  is 

0.0  =  qx  +  +  33>  -  Tx  ...  Tn  -  Tx 

Rb  Rn 

The  differential  or  algebraic  equations  governing  heat 
flow  at  each  node  are  listed  below. 

no.  i,  da) 


^*3t«3!!2Sr 
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dT2  T1  -  T2  .  T3  -  T2 
Node  2,  02^ - ^3 -  +  - 


Node  3,  C3^  =  22^  +  3*^ 


d T 


R  5 


Node  4 1  04  m  .  23^  +  25^ 


/Node  5.  05  ^  =  11^  +  215.^ 


dT6  Ta  -  T 6  .  T?  -  T6 
Node  6,  C6  ^  -  -  -R3o  -  +  YRif  - 


Koae7, 


(19) 

(20) 

(21) 

(22) 

(23) 

(24) 


....  „  „n  dT8  T7  -  T8  .  T9  -  T8  .  T16  -  T8  ,or\ 
Node  8,  C8  ^  -  -'se  +  rT3~  +  R9 U5) 


(26) 


dTIO  T9  -  T10  .  T15  -  T10 
Node  10,  CIO  +  %5 - 


ee  „ee  dTl 1  _  *  -  Til  .  T1 5  -  Til 

Node  11,  Cll  -  as  =  — ktt -  +  — nr™ - 


“ar  rt^ 


. . dT12  Aa  -  T12  .  T15  -  T12 

Node  12,  C12  sx  -  — ■  rsivr —  +  — “‘str - 


~W  =  “TSIT 


W 

L=J 

r3T 


_  =„  dTl 3  T17  -  T13  /  T14  -  T13 

Node  13,  C13  -ff - R20  + - WT^ 


(27) 

(28) 

(29) 

(30) 


Node  14i  014  ^  -  QM  +  215^114  +  11^  („, 


X** 


'wacswg gpj 
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Node  15 « 


0.0  =  QT  +  Q  +  Tl~-  Tti  +  TX6  -  T15 
u,u  WL  +  WAUX  R7  +  R8 

.  T14  -_T15  ,  T10  -  T15  .  Til  -  T15 
+  m  +  R15  Rl? 

♦  ^  <*> 


Noae  16,  0.0  = 

^  T8  -  Tl6 


R9 


T 


**•«■  015  ^  + 


(33) 

(34) 


The  Computer  Program 

The  computer  program  solves  the  differential  and 
algebraic  equations  numerically,  using  an  Euler  method 
predictor,  and  a  trapezoid  rule  corrector.  The  program¬ 
mer  selects  the  time  step  increment  and  the  total  number 
of  hours  that  are  to  be  simulated.  The  program  will  set  the 
time  of  day  to  midnight  and  then  calculate  the  temperature 
of  each  node  at  each  increment  of  time.  Print-out  includes 
nodal  temperatures,  as  well  as  the  computed  values  of  out¬ 
side  temperature,  solar  direct  gain  and  indirect  gain  heat 
inputs,  and  any  heat  added  by  the  furnace  and  internal 
equipment.  Even  though  this  study  does  not  include  heat 
from  internal  equipment,  the  provision  is  in  the  program. 

One  line  is  printed  for  each  time  increment.  A  listing  of 
the  program  along  with  an  explanation  of  the  variables  and 
information  on  required  data  input  comprizes  Appendix  A. 
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IV  Model  Verification 


At  the  Lob  Alamos  National  Laboratory  in  Lob  Alamos, 

New  Mexico,  McFarland  (Ref  8)  oversaw  a  project  to  collect 

data  on  energy  flow  in  buildings  passively  heated  by  the 

sun.  Fourteen  separate  400  ft^  cells  were  constructed, 

each  with  a  different  design  for  solar  gain.  Every  cell 

o 

had  a  collector  area  of  23.4  ft  and  faced  directly  south. 

The  cells  were  instrumented  with  thermocouples  at  important 
locations  and  equipped  with  a  600V  heater  (6  100W  light 
bulbs)  to  maintain  a  minimum  interior  temperature  of  65  F. 
Temperatures  in  the  cells,  and  the  energy  used  by  the 
heater  in  each  cell,  as  well  as  ambient  temperature,  solar 
insolation  and  wind  conditions  were  recorded  each  hour  for 
the  winter  of  1980-81. 

Test  Cell  Simulated 

Test  cell  #1  was  of  2X4  frame  construction  with  fiber¬ 
glass  and  polystyrene  insulation.  Plywood  covered  the  floor 
and  the  outside  walls.  The  cell  was  equipped  with  a  concrete, 
unvented  Trombe  wall  behind  a  double  glazed  aperture.  The 
Trombe  wall  was  15-5/8  in.  thick  with  its  south  facing  sur¬ 
face  painted  flat  black.  Forced  ventilation  maintained  an 
infiltration  rate  at  a  steady  3  complete  air  changes  per 
hour.  The  cell  was  constructed  to  have  a  Building  Load 
Coefficient  (BLC)  of  26  B/hr-F,  that  is,  to  have  a  heat 
loss  of  26  E/hr  for  every  degree  temperature  difference 
between  the  inside  and  outside.  Cell  #1  had  thermocouples 
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placed  in  the  center  of  the  interior  and  also  at  three 
locations  on  the  Trombe  walli  the  south  surface  at  mid- 
height*  the  center,  and  the  north  surface  at  mid-height. 

Cell  #1  was  modeled  so  that  these  locations  corresponded  to 
nodes  15,  1,  3#  and  5*  respectively.  Conditions  on  14 
December  1980  were  simulated.  Details  of  construction, 
along  with  calculated  input  information  for  the  computer 
program  is  given  in  Appendix  B. 

Computer  Program  Modifications 

The  computer  program  calculates  hourly  insolation 
rates  on  a  vertical  surface  given  monthly-average  daily 
total  radiation  on  a  horizontal  surface.  In  addition,  the 
program  estimates  outside  temperatures  from  monthly  averages. 
However,  data  collected  at  Los  Alamos  included  hourly  total 
radiation  on  a  vertical  surface  I^v  and  on  a  horizontal 
surface  I,  and  the  outside  temperature.  The  computer  pro¬ 
gram  was  modified  to  read  in  hourly  insolation  and  ambient 
temperature  at  each  simulation  time  step  rather  than  esti¬ 
mate  hourly  values  from  monthly  data. 

On  a  clear  day,  according  to  Duffie  and  Beckman,  the 
amount  of  diffuse  hourly  sky  radiation  1^  is  given  by 

Id  =  0.177  I  (35) 

The  hourly  beam  radiation  on  a  vertical  surface  I^y  is  the 
total  hourly  radiation  I^y  less  all  diffuse  sky  and  ground 
radiation. 

Jbv  "  *tv  -  *d  /  2  ’  fi  1  *  2 
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(36) 


The  heat  gain  at  the  outside  surface  pt  the  Trombe  wall 
is  then  given  by 

=  ^W^bv^^b  +  +  Pg  It^^r^3  /  (37) 

Results 

The  outside  temperature  measured  in  Los  Alamos,  along 
with  the  measured  interior  temperature  and  the  interior 
temperature  calculated  by  the  simulation  are  shown  in  Fig 
11.  The  actual  and  simulated  furnace  output  required  to 
maintain  the  interior  temperature  is  graphed  in  Fig  12. 

While  the  interior  temperature  was  closely  maintained  by 
the  computer  program,  difficulty  was  encountered  dupli¬ 
cating  the  furnace  behavior.  The  size  of  the  simulation 
time  step  was  0.25  hr.  When  conditions  required  that  the 
furnace  be  turned  on  in  the  simulation,  the  furnace  ran 
for  a  minimum  of  15  minutes.  The  furnace  output  rate  was 
2047  B/hr  (600  W),  so  512  B  was  input  during  each  time 
step  that  the  furnace  was  running.  This  was  far  in  excess 
of  requirements  and  the  simulated  interior  temperature  rose 
too  high  during  the  quarter  hour  that  the  furnace  was 
running.  In  the  test  cell  at  Los  Alamos,  the  heater 
cycled  as  necessary  to  maintain  the  interior  temperature, 
probably  switching  on  and  off  several  times  every  15  minutes. 
As  the  furnace  running  time  in  the  simulation  could  not  be 
decreased  without  lengthening  the  computer  processing  time 
enormously,  the  furnace  output  rate  was  decreased  by  trial 
and  error  to  900  B/hr,  presuming  that  the  furnace  would 
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Timet  hr 

Furnace  Output  Required  to  Maintain  the  Interior  Temperature 


cycle  several  times  to  generate  900  B  in  1.0  hr*  or  225  B 
in  0.25  hr.  Still,  with  only  4  time  steps  in  each  hour, 
the  simulation  had  only  4  opportunities  to  run  the  furnace 
each  hour.  Hence,  over  the  entire  hour,  only  5  different 
values  for  simulated  furnace  output  were  possible »  0,  225, 

450,  675*  and  900  B.  Despite  the  crude  approximation,  the 
total  auxiliary  energy  required  for  the  24  hr  simulation 
was  15075  B,  within  4  percent  of  the  value  measured  in 
Los  Alamos  1  14543  B. 

Temperatures  for  the  Trombe  wall  are  recorded  in  Pigs 
13,  14,  and  15.  The  simulated  thermal  behavior  of  the  wall 
was  very  similar  to  the  instrumented  wall  in  Los  Alamos. 
Temperatures  rose  and  fell  at  the  same  times  at  the  same 
locations  in  the  wall.  Simulated  temperatures,  however, 
fell  slightly  faster  during  the  early  morning  and  did  not 
rise  as  quickly  under  the  influence  of  solar  heating.  The 
maximum  difference  between  measured  and  simulated  tempera¬ 
tures  was  at  node  1,  the  south  surface  of  the  wall,  at  peak 
temperature.  The  spread  was  14  F.  The  overall  measured 
temperature  swing  was  from  63  P  to  138  F,  or  75  F.  The  sim¬ 
ulated  swing  was  from  61  F  to  124  F,  or  within  16  percent  of 
measured  values.  After  24  hours,  the  differences  at  the 
three  locations  were  4  F  at  the  south  surface  and  center, 
and  5  F  at  the  north  surface  of  the  wall. 

Fart  of  the  temperature  differences  at  node  1  can  be 
explained  by  the  simplicity  of  the  model.  The  entire  mass 
of  the  Trombe  wall  is  lumped  into  5  capacitors.  Node  1  con¬ 
tains  1/8  of  the  total  mass.  Henoe,  the  average  temperature 
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Temperature  of  the  South  Surface  of  the  Trombe  Wall 


of  the  outside  1-61/64  in.  thickness  of  the  wail  was  reported 
for  node  1,  rather  than  the  surface  temperature.  This  does 
not  explain  all  the  difference  between  simulated  and  measur¬ 
ed  temperatures,  for,  after  the  sun  went  down  the  surface 
temperature  would  drop  below  the  average  temperature  of  the 
first  1/8  thickness »  this  cross  over  did  not  take  place. 

Other  possible  sources  of  error  exist  with  the  computer 
program  and  with  the  data  taken  at  Los  Alamos.  Much  of  the 
input  information  for  the  computer  program  is  estimated.  For 
example,  R  values  for  insulation  are  nominal,  rather  than  ex¬ 
act,  and  are  affected  by  humidity.  The  resistance  due  to 
natural  convection  air  films  can  also  only  be  estimated.  The 
first  order  equation  solver  used  by  the  program  can  introduce 
errors  at  each  time  step.  In  Los  Alamos,  thermocouple  output 
can  vary.  The  pyranometers  used  to  measure  solar  radiation 
can  have  calibration  shifts.  Even  if  a  high  quality  pyrano- 
meter  is  frequently  calibrated,  errors  in  long  term  readings 
can  be  well  over  5  percent.  McFarland  did  not  discusB  pos¬ 
sible  errors  in  pyranometer  or  thermocouple  measurements. 

Refer  to  Figs  5»  13 »  and  15.  If  the  capacitance  and 
internal  resistance  of  the  Trombe  wall  has  been  accurately 
calculated,  the  slightly  faster  temperature  drop  at  nodes 
1  and  5  before  the  sun  came  up  could  be  caused  by  under¬ 
estimating  the  values  of  resistor  R2,  the  resistance  to  heat 
flow  from  the  south  surface  of  the  masonry  through  the  double 
glazing  to  the  outside,  and  resistor  R7»  the  resistance  due 
to  the  natural  convection  air  film  on  the  inside  of  the 
Trombe  wall.  To  illustrate  the  sensitivity  of  the  computer 
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program  to  the  errors  discussed  above .  the  simulation  was  re¬ 
run  with  the  values  of  resistors  R2  and  R7  increased  10  per¬ 
cent  .  and  the  value  of  hourly  solar  radiation  increased  also 
by  10  percent.  The  results  are  shown  in  Figs  16,  17.  and  18. 
There  was  no  change  in  the  furnace  output,  and  insignificant 
change  in  the  interior  temperature .  so  these  were  not  re¬ 
graphed.  The  increase  in  resistors  R2  and  R7  did  not  signi¬ 
ficantly  alter  the  behavior  of  nodes  1  and  5  before  the  sun 
came  up.  although  simulated  temperatures  approximated  measur¬ 
ed  temperatures  much  more  closely.  This  indicates  that  the 
computer  model  is  not  very  sensitive  to  errors  in  the  esti¬ 
mation  of  these  resistors.  The  marked  change  in  the  dif¬ 
ference  between  measured  and  simulated  temperatures  during 
the  period  of  solar  gain  indicates  that  the  model  is  sensi¬ 
tive  to  errors  in  solar  input.  The  maximum  difference  be¬ 
tween  simulated  and  measured  temperatures  at  node  1.  the 
south  surface  of  the  Trombe  wall,  dropped  from  14  F  to  6  F. 
The  simulated  swing  was  within  7  percent  of  the  measured 
swing.  Temperature  differences  at  the  end  of  the  simulation 
were  zero  at  the  south  surface  of  the  wall,  and  1  F  at  the 
center  and  north  surface. 

In  summary,  the  computer  model  does  a  good  job  at  pre¬ 
dicting  the  furnace  energy  required  for  a  24  hour  period, 
but  cannot  follow  hourly  trends.  In  order  to  predict  tem¬ 
perature  trends  accurately,  care  must  be  taken  in  estimating 


o 


thermal  parameters  and  inputs. 
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Pig.  18  North  Surface  of  the  Trombe  Wall  for  the  Modified  Simulation 


V  Example  Analysis 


To  Illustrate  the  usefulness  of  the  computer  program 

In  analysing  thermal  performance  of  a  passive  solar  design, 

consider  the  following  hypothetical  situation.  An  agency 

2 

is  considering  construction  of  a  1500  ft  structure  in 
Columbus,  Ohio  to  be  heated  for  office  activities.  Insula¬ 
tion  for  the  building  is  to  be  in  accordance  with  guide¬ 
lines  set  forth  in  Mazria.  The  building  is  to  be  50  ft  by 
30  ft  with  an  8  ft  ceiling.  One  $0  ft  wall  will  face 
directly  south.  R19  insulation  will  be  installed  in  the 
walls,  and  R34  in  the  ceiling.  Perimeter  insulation,  val¬ 
ued  at  RIO  will  protect  the  6  in.  thick  concrete  Blab  floor 
by  extending  2  ft  below  ground  level.  Infiltration  will  be 
kept  to  3/4  of  an  air  change  each  hour.  Window  area  will 
be  150  ft2,  distributed  evenly  around  the  walls.  All  open¬ 
ings  will  be  double  glazed.  Details  of  construction  and 
computer  input  information  is  presented  in  Appendix  C. 

To  obtain  the  annual  thermal  performance  of  this  con¬ 
ventional  structure,  a  simulation  was  run  for  an  average 
day  of  each  month.  The  total  furnace  heat  output  for  each 
day  was  multiplied  by  the  number  of  days  in  the  month,  then 
the  products  summed.  The  results,  along  with  appropriate 
heating  bills  in  Columbus,  Ohio  in  1980  for  a  gas  furnace 
(78  percent  efficient)  and  for  an  electric  furnace  (100  per¬ 
cent  efficient)  are  shown  in  Table  IV. 

As  an  alternative  to  conventional  construction,  the 
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Table  IV  Monthly  Furnace  Output  for  Conventional  Structure 


Month 

Output,  1000  B 

Gas  Bill ,  $ 

Electric 

Jan 

9610 

42.10 

130.76 

Feb 

7980 

34.96 

108.58 

Mar 

6510 

28.52 

88.58 

Apr 

3750 

16.43 

51.03 

May 

1240 

5.43 

16.87 

Jun 

0 

0 

0 

Jul 

0 

0 

0 

Aug 

0 

0 

0 

Sep 

450 

1.97 

6.12 

Oct 

1860 

8.15 

25.31 

Nov 

5550 

24.31 

75.52 

Dec 

8835 

38.70 

120.22 

Total 

45785 

200.57 

622.99 

agency  is  considering  implementing  a  rule  of  thumb  devel¬ 
oped  by  Balcomb  for  Columbus,  Ohio  which  recommends  instal¬ 
ling  a  passive  solar  collector  in  the  south  wall  with  an 
aperture  area  no  greater  than  28  percent  of  the  floor 
area.  According  to  the  rule,  larger  apertures  or  poorly 
designed  systems  with  this  large  an  aperture  would  be  prone 

to  overheating  on  a  clear  sunny  winter  day.  The  collector 

2 

area  for  a  1500  ft  floor  area  should  be  no  greater  than 
420  ft  .  The  design  under  consideration  is  a  structure 
similar  to  the  conventional  one,  but  with  the  entire  50  ft 
by  8  ft  south  wall  built  as  a  collector.  The  collector 


will  be  triple  glazed  and  will  be  partly  unvented  Trombe 

wall  and  partly  direct  gain.  Non-south  window  area  will  be 

-■  o 

75  ft  ,  double  glazed,  and  distributed  evenly  about  the  non¬ 
south  walls.  R9  night  insulation  will  be  available  for  in¬ 
stallation  on  the  outside  of  the  triple  glazing. 

Three  questions  arise  concerning  the  thermal  perform¬ 
ance  of  the  design.  How  should  the  collector  be  split  up 
between  direct  and  indirect  gain  systems  to  save  the  most 
on  the  heating  bills?  Should  night  insulation  be  used? 
Finally,  if  subjected  to  a  series  of  consistently  sunny  days, 
will  the  design  overheat?  To  find  answers  to  these  questions, 
a  simulation  was  run  for  each  month  of  the  year  for  a  col¬ 
lector  that  was  10  percent  Trombe  wall,  90  percent  direct 
gain.  Four  other  simulations  were  run  with  indirect  to  di¬ 
rect  gain  ratios  of  30/70,  50/50,  70/30,  and  90/10.  Results 
of  the  50/50  simulation  are  shown  in  Table  V.  For  compari¬ 
son,  the  empirical  SSF  method  was  applied  to  the  50/50 
simulation.  Results  are  in  the  same  table.  Annual  savings 
for  all  5  designs,  with  and  without  night  insulation  instal¬ 
led  when  the  sun  went  down  and  removed  when  the  sun  came  up 
are  graphed  in  Fig  19. 

Before  a  direct  comparison  of  the  Solar  Savings  Frac¬ 
tion  as  calculated  from  the  simulation  and  developed  by  the 
SSF  method  is  made,  the  question  of  overheating  should  be 
discussed.  The  simulations  run  to  obtain  the  results  in 
Table  V  and  Fig  19  were  for  average  days  in  each  month. 

These  days  were  neither  completely  sunny  nor  cloudy.  If 
the  design  structure  were  to  be  subjected  to  clear,  sunny 


Table  V  Furnace  Output  and  Fuel  Savings  for  the 
50/50  Passive  Solar  Design 

Without  Night  Insulation  With  Night  Insulation 

Simulation  SSF  Simulation  SSF 

Method  Method 


Month  Output 
1000  B 

Savings 

Savings 

Output 
1000  B 

Savings 

Savings 

Jan 

6975 

27 

-6 

5735 

40 

29 

Feb 

5460 

32 

0 

3780 

53 

36 

Mar 

3565 

45 

6 

2480 

62 

47 

Apr 

3750 

0 

29 

3300 

12 

67 

May 

930 

25 

66 

465 

62 

90 

Jun 

0 

- 

- 

0 

- 

- 

Jul 

0 

- 

- 

0 

- 

- 

Aug 

0 

- 

- 

0 

- 

- 

Sep 

0 

100 

100 

0 

100 

100 

Oct 

0 

100 

62 

0 

100 

88 

Nov 

3150 

43 

8 

1950 

65 

49 

Dec 

6665 

25 

-10 

5115 

42 

28 

Total 

30495 

33 

8 

22825 

50 

44 

days,  the  furnace  output  for  those  days  would  be  less  than 
for  an  average  day.  If  the  conditions  were  too  sunny,  the 
interior  temperature  might  climb  uncomfortably  high. 

As  an  example,  a  simulation  was  rerun  for  each  design, 
with  and  without  night  insulation,  for  a  24  hour  period 
following  5  consistently  sunny  days  in  January.  The  monthly 
average  dearness  index  for  January  in  Columbus ,  Ohio  is 
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0.356.  According  to  Liu  and  Jordan  (Ref  9)»  on  a  clear  day 
the  daily  total  radiation  can  be  2.4  times  greater  than  the , 
monthly-average  daily  total  radiation.  Therefore,  1102.6 
B/ft  -day,  rather  than  the  monthly-average  of  459.4  B/ft  -day 
was  input  in  each  simulation  for  the  daily  total  radiation 
H.  Figures  20  and  21  contain  plots  of  interior  temperature 
for  the  5  designs.  Table  VI  lists,  for  each  design,  the 
extreme  interior  temperature  swing,  and  for  peak  conditions, 
the  interior  air  temperature  T^,  the  mean  radiant  tempera- 
ture  T  for  an  object  in  the  center  of  the  interior,  and 
the  equivalent  uniform  temperature  T  ,  at  50  percent  rela- 
tive  humidity. 

For  office  activity  while  wearing  a  typical  business 
suit  in  50  percent  relative  humidity,  the  comfortable  tem¬ 
perature  range  from  Tables  II  and  III  is  69.5  to  80.5  F. 
Subtracting  7  degrees  from  this  range  to  account  for  in¬ 
ternally  generated  energy,  the  design  comfort  range  becomes 
62.5  to  73.5  F.  The  extreme  allowable  temperature  swing  is 
11  F.  For  simulations  of  designs  without  night  insulation, 
only  the  70/30  case  and  the  90/10  case  keep  the  temperatures 
within  the  comfort  range.  None  of  the  cases  where  night  in¬ 
sulation  is  used  have  a  peak  temperature  within  the  comfort 
range.  The  primary  conclusion  that  can  be  drawn  from  this 
information  is  that,  for  most  of  the  designs,  excess  energy 
will  have  to  be  vented  to  prevent  overheating  on  a  sunny  day. 
The  energy  dump  will  result  in  a  decrease  in  the  overall 
monthly  solar  savings. 
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Table  VI 

The  Peak 

Interior  Temperature  on 

a  Sunny 

Day 

Indirect  to 
Direct  Gain 
Area  Ratio 

Night  Tempera- 

Insulation  ture 

Swing,  P 

Tit  P 

T  »  P 
mr 

Teu’  F 

10/90 

12.8 

84.2 

82.4 

83.4 

10/90 

X 

11.2 

90.4 

88.4 

89.5 

30/70 

10.2 

80.5 

79.4 

80.0 

30/70 

X 

9.0 

87.4 

85.4 

86.5 

50/50 

7.5 

76.2 

75.7 

76.0 

50/50 

X 

6.8 

83.9 

82.7 

83.4 

70/30 

4.8 

71.4 

71.5 

71.4 

70/30 

X 

4.4 

79.6 

79.3 

79.5 

90/10 

1.9 

66.9 

67.6 

67.2 

90/10 

X 

2.1 

74.5 

75.3 

74.9 

Nonetheless,  there  is  not  a  strong  correlation  between 
the  savings  computed  by  the  simulation  and  by  the  SSF  method, 
despite  the  capability  of  the  program  to  closely  depict  the 
performance  of  the  Los  Alamos  test  cell.  Using  the  simula¬ 
tion  values  as  a  base,  the  annual  savings  predicted  by  the 
empirical  method  differed  by  as  much  as  76  percent. 
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VI  Conclusions  and  Recommendations 


The  objective  of  this  study  was  to  create  a  small 
computer  simulation  program  that  is  quick  and  easy  to  use* 
yet  can  predict  with  reasonable  accuracy  a  temperature  swing 
and  energy  savings.  It  appears  that  a  useful  tool  has  been 
devised.  Favorable  correlation  with  test  cell  experiments 
in  Los  Alamos,  New  Mexico  was  achieved.  There  is  weaker 
correlation  with  results  derived  using  the  empirical  SSF 
method.  As  pointed  out  in  Ref  3*  the  D0E2  simulation  pro¬ 
gram  also  produces  results  different  from  the  SSF  method. 

This  warrants  further  investigation.  There  are  data  for  13 
more  test  cells  in  the  Los  Alamos  study,  and  data  for  many 
more  days  than  the  day  simulated  for  test  cell  #1.  It  is 
strongly  recommended  that  this  data  be  used  for  further 
verification  studies.  Then,  the  program  should  be  used  as 
a  tool  for  a  real,  perhaps  local,  design  project. 

While  the  program  is  simple  and  requires  little  computer 
time,  it  is  not  the  best  possible  model.  Additional  features, 
such  as  the  ability  to  "open  a  window"  when  the  interior 
temperature  gets  too  high,  the  ability  to  investigate  solar 
cooling,  and  the  ability  to  specify  overhangs  for  shading 
would  make  the  program  much  more  useful.  It  is  recommended 
that  these  features  be  added  to  the  program. 

The  transient  simulation  program  was  written  for  use 
on  the  new  micro-computers.  The  final  recommendation  is 
that  a  project  to  adapt  the  program  to  such  a  piece  of  hard¬ 
ware  be  undertaken. 
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Appendix  A 

This  appendix  contains  the  computer  program  used  to 
simulate  heat  flow  in  a  structure  heated  passively  by  the 
sun.  The  program  is  written  in  PASCAL.  Following  the 
source  code  listing  is  an  explanation  of  the  significant 
variables  used  in  the  program.  Finally,  a  list  of  input 
information  in  required  order  is  presented. 
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Variables  used  Primarily  by  the  Main  Pro, 


T(l)  through  T(l?)i  the  temperatures  at  nodes  1  through 
17  at  the  current  time. 

RUNTIME 1  the  number  of  hours  that  the  program  is  to 
simulate.  The  minimum  value  of  this  input  variable  is 
264.0,  for  the  program  runs  internally  for  240.0  simulated 
hours  to  reach  cyclic  steady  state. 

FURNACERATEi  the  heat  output  per  hour  of  the  furnace. 

LOWSETt  the  temperature  below  which  the  furnace  kicks  on 

HISETi  the  temperature  above  which  the  furnace  turns  off 

TAVEi  the  monthly- average  temperature. 

TSWINGi  half  the  difference  between  the  monthly-average 
high  and  the  monthly-average  low  temperature. 

TIME 1  the  current  time. 

NEXTTIME 1  the  time  at  one  time  step  into  the  future. 

OPENTIME 1  the  time  of  day  that  night  insulation  is  to 
be  removed. 

CLOSETIMEi  the  time  of  day  that  night  insulation  is  to 
be  installed. 

MONTH t  an  integer  flag  to  identify  the  month  of  the  year 
This  input  variable  will  be  assigned  a  value  from  1  to  12. 

RINi  the  value  of  resistor  R1  when  night  insulation  is 
not  installed. 

NITEFACTORt  the  value  by  which  R1N  is  multiplied  when 
night  insulation  is  in  place. 

R9Ni  the  value  of  resistor  R9  during  the  day. 


SMNITEFACTORi  the  value  by  which  R9  is  multiplied  when 
the  system  is  operating  in  night  configuration. 

Variables  used  Primarily  by  the  Differential  Equation 

Solver. 

TP(1)  to  TP(l?)i  the  predicted  temperatures  at  the  nodes. 

TC(1)  to  TC(17)»  the  corrected  temperatures  at  the  nodes. 

TI(1)  to  TI(17)»  a  dummy  variable  used  to  hold  calculated 
nodal  temperatures. 

DT(1)  to  DT(l7)i  the  temperature  gradients  at  the  nodes. 

DT0(1)  to  DT0(17)i  a  dummy  variable  used  to  hold  calcu¬ 
lated  values  of  temperature  gradients. 

TAUi  the  time  of  day. 

DTAUt  the  time  step  size  in  hours.  In  R-C  circuit 
analysis,  the  time  step  size  is  determined  from  the  relation! 

DTAU  «  Minimum  (Cx/£(1/Ry)) 

where  Cx  is  the  value  of  each  capacitor  and  Ry  is  the  value 
of  each  resistor  directly  connected  to  capacitor  Cx. 

RELERRt  the  relative  error  between  the  pi . Ucted  and  the 
corrected  temperature  at  each  node. 

EPSILON i  the  maximum  allowable  relative  error.  A  value 
of  0.001  is  recommended. 

TAi  the  ambient  temperature  at  the  current  time  of  day. 

TANEXTi  the  ambient  temperature  at  one  time  step  in  the 
furure. 

TAI i  a  working  variable  for  ambient  temperature. 

QAOXNEXTi  the  furnace  input  for  node  15  at  one  time 
step  in  the  future. 
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QLOADNEXT t  the  internally  generated  heat  input  for 
node  15  at  one  time  step  in  the  future. 

R1  to  R25»  the  values  of  resistors  R1  to  R25. 

Cl  to  C15«  the  values  of  capacitors  Cl  to  C15. 

Variables  used  in  the  Calculation  of  Solar  Gain. 

CNi  the  total  number  of  glazings  on  the  collector. 

NCt  the  index  of  refraction  of  glass. 

KAPPA «  the  extinction  coefficient  of  glass. 

THICKGLASs  the  thickness  of  one  of  the  glazings  on  the 
collector. 

Al,  A2,  and  A31  collector  loss  ratio  coefficients. 

ALFA:  the  effective  absorptance  of  the  interior  surfaces. 
ALFABs  the  beam  absorptance  of  the  surface  of  the  Trombe 

wall. 

ALFABN  t  the  normal  beam  absorptance  of  the  surface  of 
the  Trombe  wall. 

ALFADi  the  diffuse  absorptance  of  the  surface  of  the 
Trombe  wall. 

RHDGs  the  reflectance  of  the  ground  surrounding  the 
collector. 

RHORt  the  effective  reflectance  of  the  interior  surfaces. 
RHD1,  RH02 1  the  reflectance  of  the  collector  glazings. 
TAUBNt  the  normal  beam  transmittance  of  the  collector 
glazings. 

TAUBi  the  beam  transmittance  of  the  collector  glazings. 

TAUDi  the  diffuse  transmittance  of  the  collector  glazings 

TAB 1  the  effective  transmit tance-absorptance  product  for 

beam  radiation  inoident  on  the  Trombe  wall. 
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TAD«  the  effective  transmittance-absorptance  product 
for  diffuse  radiation  incident  on  the  Trombe  wall. 

TAUAi  the  transmittance  due  to  absorptance  of  all  of 
the  collector  glazings. 

TAUAI i  the  transmittance  due  to  absorptance  of  one  of 
the  collector  glazings. 

TAUli  the  total  transmittance  of  one  of  the  collector 
glazings. 

DGAREAt  the  area  of  the  collector  dedicated  to  direct 

gain. 

TWAREAi  the  area  of  the  collector  dedicated  to  indirect 

gain. 

rOOMAREA  t  the  surface  area  of  the  interior,  less  glazed 

area. 

QTWi  the  indirect  gain  heat  input  at  node  1  at  the  current 

time. 

QTWNEXT i  the  indirect  gain  heat  input  at  node  1  at  one 
time  step  in  the  future. 

Q1MI,  TWt  working  variables  for  indirect  gain  heat  inputs. 
QDGt  the  direct  gain  heat  input  at  node  14  at  the  current 

time. 

QDGNEXTt  the  direct  gain  heat  input  at  node  14  at  one 
time  step  in  the  future. 

QD6I,  DGi  working  variables  for  direct  gain  heat  inputs. 
OURANGL i  the  hour  angle. 

SSOURANGLi  the  sunset  hour  angle. 

THETA i  the  angle  of  incidence  of  beam  radiation. 


SINTHETA «  sin( THETA ) . 

.  COSTHETAt  cos (THETA). 

THETACi  the  angle  of  beam  radiation  within  the  glazing 
after  refraction  at  the  air-glass  interface. 

SINTHETACi  sin(THETAC) , 

COZENITHANGLi  the  cosine  of  the  angle  of  incidence  of 
beam  radiation  on  a  horizontal  surface. 

DAYOFYEARi  Julian  date. 

AZMTHi  azimuth  in  degrees. 

LATi  latitude  in  degrees. 

DECLNATN :  declination. 

RBt  the  ratio  of  beam  radiation  on  a  tilted  surface  to 
beam  radiation  on  a  horizontal  surface. 

HT:  the  daily,  total  radiation  on  a  horizontal  surface. 
HDi  the  daily  diffuse  radiation  on  a  horizontal  surface. 
KTi  the  clearness  index. 

IT:  the  hourly  total  radiation. 

ID«  the  hourly  diffuse  radiation. 

IB i  the  hourly  beam  radiation. 

Required  Input  Information 

Values  for  the  following  variables  are  required  for 
the  input  data  file.  Information  must  be  in  the  order  listed 
below,  b’vtt  muy  be  grouped  in  any  way.  Asterisked  values  are 
integers.  All  others  are  real  numbers. 

Parameters  of  the  Structure  and  Location. 

R1N,  NITEFACTOR,  R2,  R3.  R4,  R5t 
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R6,  R7t  R8,  R9N,  SMWITEFACTOR,  RIO, 

Rll,  R15,  R16,  R17i  R18,  R19, 

R20,  R22,  R23,  R24,  R25. 

Cl,  C2,  C3,  C4,  C5t  C7,  Cll,  C12,  C13, 

LAT,  AZMTH,  CN*,  Al,  A2,  A3, 

THICKGLAS ,  NC,  KAPPA,  ALFA BN,  RHOR,  RHOG, 

DGAREA,  3WAREA,  ROOMAREA ,  FURNA CERATE,  LOWSET,  HISET, 
Particulars  about  the  Current  Run. 

MONTH*,  RUNTIME,  DTAU,  EPSILON,  DAYOFYEAR*, 

OPENTIME,  CLOSETIME,  HT,  KT,  TAVE,  TSWING, 

Initial  Conditions. 

TA ,  QLOADNEXT,  T(l),  T(2),  T(3).  T(4), 

T(5).  T(6),  T(7),  T(8),  T(9) f  T(10),  T(ll), 

T(  12 ) ,  T(13),  T(14),  T(15),  T(16,  T(17). 

If  the  program  is  modified  to  utilize  internally  generated 
heat  Ql,  the  hourly  heat  input  rate  must  be  included  in 
the  input  file  for  each  time  step  in  the  simulation! 
QLOADNEXT,  QLOADNEXT,  .... 
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Appendix  B 


Details  of  the  construction  of  test  cell  #1  in  Los 
Alamos,  New  Mexico,  along  with  calculated  input  information 
for  the  modified  computer  program  is  presented  in  this 
appendix.  Figure  4  and  Appendix  A  will  he  useful  in  follow¬ 
ing  the  development  here.  Unless  otherwise  indicated,  mater¬ 
ial  properties  are  from  Mazria  or  the  Carrier  System  Design 
Manual  (Ref  10). 

The  test  cell  measured  86  in.  (north-south)  by  62  in. 
(east-west)  by  120  in.  The  west  wall  of  the  cell  was  a  com¬ 
mon  wall  between  cell  #1  and  cell  #2.  Because  the  wall  was 
well  insulated  and  because  temperatures  on  both  sides  of 
the  wall  were  maintained  at  or  near  65  F#  it  was  assumed 
that  no  heat  transfer  took  place  across  the  wall. 

The  test  cell  did  not  have  direct  gain  heating,  nor 
a  masonry  floor,  a  detached  storage  mass,  or  a  sunspace. 

For  this  reason,  nodes  6,  7,  8,  9*  10,  13»  14,  16,  and  17 
were  eliminated  from  the  circuit.  This  was  done  by  setting 
the  values  of  resistors  R8,  R9.  RIO,  Rll,  R12,  R13»  R14,  R15, 
R20,  R21,  R24,  and  R25  to  1X1015  hr-F/B,  resistor  R1  to  1X10~15 
hr-F/B,  and  capacitors  C6,  C7*  C8,  C9«  CIO,  C13,  C14,  and  C15 
to  0.01  B/F.  In  addition,  the  test  cell  had  windows  and  no 
significant  thermal  mass  in  the  conventional  walls.  For 
this  reason,  resistor  R22  was  used  to  model  resistance  to  heat 
flow  in  the  walls,  and  node  11,  with  associated  resistors  R16 
and  R17  and  capacitor  Cll  was  used  for  the  floor  of  the  cell. 
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Resistor  R2 .  The  U-value  for  the  double  glazings  and 
air-space  in  front  of  the  Trombe  wall  was  given  in  McFarland 
as  0.46  B/ft2-hr-F.  The  glazed  area  is  23.4  ft2. 

R2  ■  1  /(0.46  X  23.4) 

=  0.0929023  hr-F/B 

Resistors  R3  through  R6.  From  McFarland,  the  conduc¬ 
tivity  of  th$  concrete  in  the  Trombe  wall  was  1.0  B/ft-hr-F. 
For  a  wall  with  a  cross  sectional  area  of  29.8  ft2  and  a 
thickness  of  15-5/8  in. ,  one  fourth  the  total  internal  re¬ 
sistance  is  1 

R3  «  1  X  ( 15-5/8 )/( 12  X  4  X  29.8) 

=  0.0109235  hr-F/B 

R3  =  R4  «  R5  =  R6 

Resistor  R7.  The  natural  convection  air  film  coeffic¬ 
ient  for  a  vertical  wall  is  0.68  ft2-hr-F/B.  The  exposed 
area  of  the  Trombe  wall  is  29.8  ft2. 

R7  =  0.68/29.8 

=  0.0224161  hr-F/B 

Resistors  R16  and  R17.  The  bottom  of  the  test  cell 
sat  on  2X4  framing  and  was  exposed  to  outside  air  flow. 
Construction,  from  outside  to  inside  with  associated  re- 
sistances  in  ft  -hr-F/B  wasi 


Outside  Air  Film  0.25 

5/8  in.  Plywood  0.78125 

2X6  Frame  and  R22  Insulation  22.0 

5/8  in.  Plywood  0.78125 

Total  Resistance  23.8125 


Floor  area t  37.03  ft2 
R16  =  23.8125/37.03 
=  0.6430983  hr-F/B 
For  resistor  R17i 

Inside  Air  Film  0.92 

R17  =  0.92/37.03 
=  0.0248462 

Resistors  R18  and  R19.  Construction  of  the  ceiling  and 
roof,  from  outside  to  inside  with  associated  resistances  in 


ft2-hr-F/B  was  1 

Outside  Air  Film  0.17 

Rolled  Roofing  0.15 

90#  Felt  0.06 

5/8  in.  Plywood  0.78125 

Air  Space  0.85 

2X6  Frame  and  R22  Insulation  22.0 

Galvanized  Sheet  Metal  q.q 

Total  Resistance  24.01125 

Ceiling  Areas  37.03  ft2 
R18  «  24.01125/37.03 
■  0.6484659  hr-F/B 
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For  resistor  R19« 


Inside  Air  Film  0.61 

R19  =  0.61/37.03 

=  0.0164741  hr-F/B 

Resistor  R22.  Three  different  constructions  with  as¬ 
sociated  resistances  RA ,  RB,  and  RC  were  modeled  with  resis¬ 
tor  R22.  Construction,  from  outside  to  inside  with  associ¬ 
ated  resistances  in  ft^- hr-F/B  follow.  For  the  north  and 


east  wallt 

Outside  Air  Film  0.17 

5/8  in.  Plywood  0.78125 

2X4  Frame  and  Rll  Insulation  11.0 

1  in.  Polystyrene  5- 0 

Inside  Air  Film  0.68 

Total  Resistance  17.63125 

Area.  105.3  ft2 


RA  *  17.63125/105.3  =  0.1673853 
For  the  access  door. 


Outside  Air  Film  0.17 

5/8  in.  Plywood  0.78125 

4  in.  Polystyrene  20.0 

Inside  Air  Film  0.66 

Total  Resistance  21.63125 

RB  «  21.63125/  18.0  *  1.2017361 
For  the  south  wall. 

Outside  Air  Film  0.17 

5/8  in.  Plywood  0.78125 
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V 


22.5 

0.68 

24. 13125 


4-1/3  in.  Polystyrene 
Inside  Air  Film 
Total  Resistance 
Areat  21.87  ft2 
RC  =  24.13125/21.87  =  1.1035633 
Resistor  R22  was  the  equivalent  of  these  three  parallel 
resistors. 

R22  =  l/( 1/RA  +  1/RB  +  1/RC) 

»  0.1296592  hr-F/B 

Resistor  R23.  The  infiltration  rate  was  3  air  changes 


per  hour. 

Applying  Eq  13 1 

p% 

0.013  B/ft3-F 

V 

370.28  ft3 

n 

3  1/hr 

Q  *  0.013  X  370.28  X  3  *  14.4408  B/hr-F 
R23  =  1/14.4408 

=  0.0692481  hr-F/B 

Capacitors  Cl  through  C5.  From  McFarland,  the  total 
capacitance  of  the  Trombe  wall  was  952  B/F. 

Cl  =  952.0/8  =  119.0  B/F 
C2  =  952.0/4  =  238.0  B/F 
C3  =  C4  =  C2 
C5  *  Cl 

Capacitor  Cll.  Only  mass  inside  the  insulation  was 
effective  as  thermal  mass.  As  an  estimation,  all  of  the 
plywood  floor  and  half  of  the  2X6  frame  running  east-west 


Volume  of  Plywood 
'  Density 
Specific  Heat 
Heat  Capacity  of  Plywood 


1.93  ft3 
34.0  Xba/tt* 
0.29  B/lbm-F 
19.02  B/P 


Volume  of  half  of  Frame 
Density 
Specific  Heat 
Heat  Capacity  of  Joists 
The  value  of  capacitor  Cll  will  be 


0.97  ft3 
27.0  lbin/ft3 

14.93  B/P 

the  sum  of  the  above 


components. 


Cll  =  14.93  +  19.02 


=  33.95  B/P 

Capacitor  C12.  Significant  mass  in  the  ceiling  comes 
from  the  sheet  metal  and  half  of  the  2X6  joists. 


Volume  of  Metal 
Density 
Specific  Heat 
Heat  Capacity  of  Metal 
Heat  Capacity  of  Joists 
C12  »  5.66  +  14.93 
»  20.59  B/P 


0.0964  ft3 
489.0  lbm/ft3 
0.12  B/lbm-F 
5.66  B/F 
14.93  B/P 


LAT.  35.8  degrees. 
AZMTH.  0.0  degrees. 
£2(.  2  gift  lings. 

41.  0.15. 

44.  0.62. 
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Cl 


< 


A2.  0.0 

,  THICKGLAS.  0.1875  in. 

NC.  1.526. 

KAPPA.  0.5452  in. 

ALPABW.  0.95. 

RHOR.  This  variable was  not  used  because  the  cell  did 
not  have  direct  gain  heating.  The  value  was  set  at  0.5. 

RHOG.  0.3  estimated. 

PGA REA.  0.0  ft2. 

TWAREA.  23.4  ft2. 

ROOMAREA.  This  variable  was  not  used  because  the  test 
cell  did  not  utilize  direct  gain.  The  value  was  set  at 
100  ft2. 

FURNA CERATE.  900  B/hr. 

LOWSET.  65  P. 

HISET.  66  P. 

DTAP.  In  R-C  circuit  analysis »  the  time  step  size  is 
determined  from  the  relation! 

DTAU  5  Minimum  (Cx/£(l/Ry)) 

where  Cx  is  the  value  of  each  capacitor  and  Ry  is  the  value 
of  each  resistor  directly  connected  to  capacitor  Cx.  At 
node  12 f  resistors  R18  and  R19  Are  attached  to  capacitor  C12. 
C12/(l/R18  +  1/R19)  -  20.59/(1/0.648  ♦  1/0.0165) 

«  0.33  hr 

This  relation  was  computed  for  all  capacitors.  Rode  12 
yielded  the  smallest  time.  As  DSfcU  must  be  smaller  than 
‘tills  minimum, 
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DTAU  =  0.25  hr 

EPSILON.  0.001. 

DAYOPYEAR.  348. 

OPENTIME  and  CLOSETIME.  As  no  night  insulation  was 
specified*  these  variables  were  assigned  values  so  that 
they  did  not  interfere  with  the  program. 

OPENTIME  =  -5.0 

CLOSETIME  =30.0 

KT.  0.75  estimated. 

TA.  28  P. 

TiU.  67.0  P. 

T(2).  69.0  P. 

T(3).  71.0  P. 

T(4).  70.5  F. 

T(5).  70.0  P. 

T(6)  to  R(10)  and  T(13).  T(14)  and  T(17).  These  cor¬ 
respond  to  isolated  nodes  and  were  set  at  99*0  P. 

T(ll).  68.0  P. 

T(  12 ) .  68.0  P. 

T(15).  68.0  P. 

T(16).  28.0  P. 

TA.  IT.  and  ITy.  The  modified  computer  program  read 
in*  at  each  time  step*  the  ambient  temperature*  and  the 
hourly  rates  of  solar  insolation  on  a  horisontal  and  on  a 
vertical  surface  in  t/tX?-hr,  These  data  were  recorded  at 


Los  Alamos  for  the  entire  1980-81  winter.  On  14  December, 


1  L5 


Time 

TA 

IT 

ITV 

01 

27 

0 

0 

02 

26 

0 

0 

03 

25 

0 

0 

04 

26 

0 

0 

05 

27 

0 

0 

06 

26 

0 

0 

07 

27 

0 

0 

08 

27 

4 

3 

09 

27 

34 

24 

10 

32 

130 

213 

11 

39 

148 

276 

12 

42 

167 

296 

13 

43 

167 

295 

14 

45 

147 

273 

15 

46 

109 

232 

16 

46 

59 

173 

17 

42 

7 

46 

18 

37 

0 

0 

19 

34 

0 

0 

20 

33 

0 

0 

21 

33 

0 

0 

22 

33 

0 

0 

23 

31 

0 

0 

24 

30 

0 

0 

ttLfegMpg.  Ifeg-JSiUdAlMLiaafl  Acoordin*  to 


cell  vu  26.0  B/hr-f  with  an  infiltration  rate  of  3/hr,  and 
a  BLC  of  16.0  B/hr-F  with  an  infiltration  rate  of 
1/hr.  Calculating  the  BLC  using  the  values  of  resistors 
estimated  above,  values  of  25. 15  and  15*53  B/hr-F  we„*e  ob¬ 
tained.  By  changing  the  value  of  resistor  R22  frost  0.1296592 
to  0.1250163  hr-F/B  and  the  value  of  resistor  R23  from 
0.0692481  to  O.O666667  hr-F/B,  the  design  BLC *8  were 
closely  matched. 


Details  of  the  construction  of  the  conventional  design 
for  Columbus ,  Ohio,  along  with  calculated  input  information 
for  the  computer  program  is  presented  in  this  appendix. 
Figure  4  and  Appendix  A  will  be  useful  in  following  this 
development. 

The  conventional  design  does  not  have  a  sunspace, 
a  Trombe  wall,  or  a  detached  storage  mass.  For  this  reason 
nodes  1,  2,  3,  4,  5,  6,  7,  8,  9,  and  10  are  eliminated  from 
the  circuit.  This  is  done  by  setting  the  values  of  re¬ 
sistors  R2,  R3,  R4,  R 5,  R6,  R?,  R9.  RIO,  Rll,  R12,  R13, 

R14,  and  R15  to  1X101-*  hr-F/B,  the  value  of  resistor  Rl  to 
1X10-1-’  hr-F/B,  and  the  values  of  capacitors  Cl,  C2,  C3.  C4 
C5,  C6,  C8,  C9,  and  CIO  to  0.01  B/F.  The  values  of  remain¬ 
ing  input  information  is  calculated  from  Mazria  or  the 
Carrier  System  Design  Manual,  as  applicable,  unless  other¬ 
wise  specified. 

Resistor  R8.  Windows  are  double  glazed  with  a  1/2  in. 

2 

air  gap  between  glazings.  With  150  ft  total  window  area 
distributed  evenly  about  the  structure,  46.875  ft  will 
be  built  into  the  south  wall. 

Transmission  Coefficient  0.55  B/ft  -hr-F 

R8  =  1/(0. 55  X  46.875) 

«  0.0387879  hr-F/B 


Resistors  R16  and  R17.  The  walls  are  constructed  from 
outside  to  in  of  wood  shingles  on  top  of  5/16  in.  plywood 
sheathing,  2X6  frame  on  24  in.  centers  with  R19  insulation, 
and  3/4  in.  wood  panel  on  the  interior.  The  wall  area  is 
the  product  of  the  wall  height  and  the  length  of  the  peri- 
meter,  less  the  window  area  or  1130  ft  .  For  resistor  RI61 


Wall  Resistance 
Less  the  Inside  Air  Film 
Total  Resistance 
R16  =  22.32/1130.0 

=  0.0197522  hr-F/B 
For  resistor  Rl?i 

The  Inside  Air  Film 


23.0 

0.68 

22.32  ft2 -hr-F/B 


0.68  ft  -hr-F/B 


R17  *  0.68/1130.0 

■  0.0006018  hr-F/B 

Resistors  R18  and  R19.  The  ceiling  and  roof  are  con¬ 
structed,  from  outside  to  in,  of  wood  shingles  on  top  of 
building  paper  and  5/16  in.  plywood  sheathing,  2X6  frame 
on  24  in.  centers  with  R34  insulation,  a  25/32  in.  wood 
sub-floor  to  which  is  glued  3/4  in.  thick  acoustical  tile. 


For  resistor  RI81 

Ceiling  Resistance 
Less  the  Inside  Air  Film 
Total  Resistance 
Floor  Area.  1500  ft2 
R18  -  39.675/1500.0 
-  0.0264501  hr-F/B 


40.285 

0.61 

39.675  ft2-hr-F/B 
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For  resistor  R19« 


_  Inside  Air  Film  0.61  ft2-hr-F/B 

R19  *  0.61/1500.0 

=  0.00040667  hr-F/B 

Resistors  R20,  R21.  R24.  and  R 25.  The  floor  is  a  6  in. 
thick  concrete  slab  overlaid  with  linoleum.  RIO  insulation 
extends  2  ft  below  ground  around  the  160  ft  perimeter.  For 
resistors  R20  and  R21t 

Resistance  of  3  in.  thick  concrete  0.24  ft2-hr-F/B 
Floor  Areat  1500  ft2 
R20  =  0.24/1500.0 
=  0.00016  hr-F/B 
R21  »  R20 
For  res is tor  R24i 

Inside  Air  Film  0.61  ft2-hr-F/B 

R24  *  0.61/1500.0 

=  0.0004066?  hr-F/B 
For  resistor  R25,  apply  Eq  12 1 

Edge  Lobs  Factor  0.17  B/ft-hr-F 

Perimeter  160.0  ft 

R25  *  1/(0. 17  X  160,0) 

*  0.0367647  hr-F/B 

Resistor  R23.  The  infiltration  rate  in  the  structure 
is  estimated  at  3/4  change  per  hr. 

0.018 

12000.0  ft3 
0.75  1/hr 


Heat  Capacity  of  Air 
Volume  of  Interior 
Air  Changes 


496.08  B/P 


18.3  ft3 
50.0  lbn/ft3 
0.31  B/lbm-F 


283.64  B/F 


R23  =  1/(0.018  X  12000.0  X  0.75) 

»  0.0617284  hr-F/B 

Capacitor  Cll.  Only  mass  inside  the  insulation  is 
effective  as  thermal  mass.  As  an  estimate ,  all  of  the 
paneling  and  half  of  the  2X6  frame  is  included. 

Volume  of  Paneling  50.3125  ft3 

Density  34.0  lbm/ft3 

Specific  Heat  0.29  B/lbm-P 

Heat  Capacity  of  Paneling 
Volume  of  Framing 
Density 
Specific  Heat 
Heat  Capacity  of  Framing 
Cll  «  496.08  +  283.64 
=  780.24  B/F 

Capacitor  C12.  All  of  the  sub- floor  and  acoustical 
tiley  and  half  of  the  frame  in  the  roof  will  be  included 
as  thermal  mass. 

Volume  of  Tile 
Density 
Specific  Heat 
Heat  Capacity  of  the  Tile 
Volume  of  Sub-Floor 
Density 
Specific  Heat 

Heat  Capacity  of  the  Sub-floor 


93.75  ft3 
18.0  ltn/ft3 
0.19  B/lbm-F 

97.66  ft3 
50.0  lbm/ft3 
0.31  B/lbm-F 


320.62  B/F 


1513.73  ** 


Volume  of  Frame 
Density 


41.6?  ft3 
50.0  lbn/ft3 


Specific  Heat  _ 

Heat  Capacity  of  Framing 
C12  =  320.62  +  1513.73  +  645.88 
*  2480.23  B/F 

Capacitors  C13.  C14.  and  CIS. 


0.31  B/lbm-F 


645.88  B/F 


Volume  of  Concrete 
Density 
Specific  Heat 
Total  Heat  Capacity 


750.0  ft3 
140.0  lbm/ft3 
0.22  B/lbm-F 


23100.0  B/F 


C13  =  11550.0  B/F 
C14  =  5775-0  B/F 
C15  =  5775.0  B/F 
LAT.  4 0.0  degrees. 

AZMTH.  0.0  degrees. 

CN.  2  glazings. 

Al.  0.14. 

A2.  0.45. 

A2.  0.75. 

THICKGLAS.  0. 125  in. 

S£.  1.526. 

KAPPA.  0.5452. 

ALPABN.  0.94. 

RHOR.  0.4. 

RHOG.  The  ground  reflectance  depends  upon  are rage  snow 
cower  for  the  month.  RHOG  was  estimated  at  0.45  for  January* 
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February,  and  December,  0.25  for  March  and  November,  and 
0.2  for  the  remaining  months. 

PGA REA.  46.875  ft2. 

TWA REA.  0.0  ft2. 

ROOMAREA.  4130  ft2. 

FURNA CERATE.  20000.0  B/hr. 

LOWSET.  65.0  F. 

HISET.  69.0  F. 

MONTH.  For  January,  1.  For  December,  12. 

RUNTIME.  To  eliminate  errors  caused  by  poor  selection 
of  initial  conditions,  240  hours  are  simulated  internally 
before  printout  begins. 

RUNTIME  «  336.0  hr. 

DTAU.  The  calculations  described  in  Appendix  B  are  carried 
out  for  each  capacitor  in  this  circuit.  The  minimum  value 
obtained  is  0.63  hr. 

DTAU  =  0.25  hr. 

EPSILON.  0.001. 

DAYOFYEAR.  Typical  days  of  each  month,  with  their 
Julian  dates,  arei  17  January  (17)»  16  February  (47),  16 
March  (75).  15  April  (105).  15  May  (135).  11  June  (162), 

17  July  (198),  16  August  (228),  15  September  (256),  15 
October  (288),  14  November  (318),  and  10  December  (344). 

OPENTIME.  -5.0. 

CLOSETIME.  30.0 

HT.  ET.  TAVE.  TSWING.  and  TA.  For  Columbus,  Ohio 
insolation  and  temperature  data  are  extracted  from  Liu 
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and  Jordan,  and  frpm  Cinquemani  (Ref  11).  Values  for  TA  are 
calculated  using  Eq  17.  Values  for  HT  are  in  B/f t2-day. 
Values  for  temperatures  are  in  degrees  F. 


Month 

HT 

KT 

TSAVE 

TSWING 

TA 

Jan 

459.4 

0.356 

28.4 

8.0 

21.5 

Feb 

676.8 

0.401 

30.3 

8.9 

22.6 

Mar 

979.6 

0.447 

39.2 

10.1 

30.5 

Apr 

1352.9 

0.470 

51.2 

11.7 

41.4 

May 

1646.9 

0.515 

61.1 

11.8 

50.9 

Jim 

1812.7 

0.561 

70.4 

11.5 

60.4 

Jul 

1754.9 

0.555 

73.6 

11.2 

63-9 

Aug 

1640.6 

0.475 

71.9 

11.8 

61.7 

Sep 

1261.6 

0.433 

65.2 

12.5 

54.7 

Oct 

945.1 

0.441 

54.2 

12.2 

43.6 

Nov 

537.8 

0.302 

41.7 

9.3 

33.6 

Dec 

387.2 

0.351 

30.7 

10.6 

21.5 

QLOADNEXTi  0.0. 

T(l)  to  T(17).  The  initial  temperatures  of  the  nodes 
are  set  at  60.0  F. 
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